Not to appear in Nonlearned J., 45. 

Preprint typeset using I^T^X style cmulatcapj v. 6/6/04 



SPITZER-IRS HIGH RESOLUTION SPECTROSCOPIC SURVEY OF THE 12 /jM SEYFERT GALAXIES: 

II. RESULTS FOR THE COMPLETE DATASET 

Silvia Tommasin 1 , Luigi Spinoglio 

Istituto di Fisica dello Spazio Interplanetario, INAF, Via Fosso del Cavaliere 100, 1-00133 Roma, Italy 



Matthew A. Malkan 

Astronomy Division, University of California, Los Angeles, CA 90095-1547, USA 

AND 

Giovanni Fazio 

Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
Not to appear in Nonlearned J., 45- 

ABSTRACT 

We present our Spitzer IRS spectroscopic survey from IO^to to 37 fim of the Seyfert galaxies of the 
12/im Galaxy Sample, collected in high resolution mode (R~ 600). The new spectra of 61 galaxies, 
together with the data we already published, gives us a total of 91 12/im Seyfert galaxies observed, out 
of 112. We discuss the mid-IR emission lines and features of the Seyfert galaxies, using an improved 
AGN classification scheme: instead of adopting the usual classes of Seyfert l's and Seyfert 2's, we use 
the spectropolarimetric data from the literature to divide the objects into categories "AGN 1" and 
"AGN 2" , where AGN l's include all broad-line objects, including the Seyfert 2's showing hidden broad 
lines in polarized light. The remaining category, AGN 2's contains only Seyferts with no detectable 
broad lines in cither direct or polarized spectroscopy. We present various mid-IR observables, such 
as ionization-sensitive and density-sensitive line ratios, the PAH 11.25/zm feature and the H2 S(l) 
rotational line equivalent widths, the (60/ito - 25/ito) spectral index and the source extendedness at 
19/xm, to characterize similarities and differences in the AGN populations, in terms of AGN dominance 
versus star formation dominance. 

We find that the mid-IR emission properties characterize all the AGN l's objects as a single family, 
with strongly AGN-dominated spectra. In contrast, the AGN 2's can be divided in two groups, the 
first one with properties similar to the AGN l's except without detected broad lines, and the second 
with properties similar to the non-Seyfert galaxies, such as LINERs or starburst galaxies. 

We computed a semianalytical model to estimate the AGN and the starburst contributions to the 
mid-IR galaxy emission at 19/im. For 59 galaxies with appropriate data, we can separate the 19/im 
emission into AGN and starburst components using the measured mid-IR spectral features. We use 
these to quantify the brightness thresholds that an AGN must meet to satisfy our classifications: AGN 
1 have an AGN contribution > 73 % and AGN 2 > 45 % of their total emission at 19/zra. 

The detection of [NeV] lines turns out to be an almost perfect signature of energy production by an 
AGN. Only 4 (~ 7.5% percent) of 55 AGN 1 and 2 (10% percent) out of 20 AGN 2 do not have [NeV] 
14.3/mi down to a flux limit of ~ 4 x 10~ 15 er<?s _1 cm -2 . We present mean spectra of the various AGN 
categories. Passing from AGN-dominated to starburst-dominated objects, the continuum steepens, 
especially at wavelengths shorter than 20fim, while the PAH feature increases in its equivalent width 
and the high ionization lines decrease. 

We estimate H2 mass and excitation temperature through the measurement of the S(l) rotational 
line of this molecule. Finally we derive the first local luminosity functions for the brighest mid-infrared 
lines and the PAH feature at 11.25/xto. No statistical difference is apparent in the space densities for 
Seyfert l's and 2's of a given line luminosity, nor for the new classes of AGN l's and 2's. We use 
the correlation between [Ne V] line and nonstellar infrared continuum luminosity to derive the global 
output of accretion-powered galactic nuclei in the local universe. 
Subject headings: Galaxies: Active - Galaxies: Starbursts - Infrared: Galaxies 



1. INTRODUCTION 

This paper contains the final results of the Spitzer 
high-resolution IRS spectroscopic survey of the sample 
of Seyfert galaxies (hereafter 12MSG) included in the 
IRAS 12/im galaxy sample (Rush, Malkan & Spinoglio 
1993, hereafter RMS). In Tommasin et al. (2008) (here- 

1 also at the Physics Department of Universita di Roma, La 
Sapicnza, Roma, Italy 



after Paper I) we have presented and analyzed the first 
30 high-resolution spectra of 29 Seyfert galaxies of this 
sample (one IRAS galaxy, Mrk 1034, was coincident 
with a pair, for which we obtained two spectra). The 
first spectroscopic observations of active galaxies with 
the Infrared Spectrometer (Houck et al. 2004, IRS) 
onboard the Spitzer Space Observatory (Werner et al. 
2004) have been collected on classical AGNs (Weedman 
et al. 2005) and ULIRGs (Armus et al. 2007). Af- 
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ter the work presented in Paper I and the referenced 
works therein, a few more studies have discussed the 
Spitzcr mid-infrared spectra of Seyfert galaxies, among 
these Deo et al. (2007), Melendez et al. (2008a) and 
Wu et al. (2009). As expected, the mid- infrared spec- 
tra of Seyfert galaxies show forbidden lines originating in 
the Narrow Line Regions, excited by the AGN ionizing 
flux. This power is thought to be produced by black-hole 
accretion, i.e., ultimately from the conversion of gravi- 
tational into radiative energy. The fine structure lines 
of [NeV] at 14.32/xm and 24.31/xm originate exclusively 
in the highly ionized gas (with an ionization potential 
of 97eV) illuminated by the AGN 2 . As discussed in Pa- 
per I, the [OIV] line at 25.88/im (ionization potential of 
55eV) is most probably excited from the AGN. In fact, 
Melendez et al. (2008b) consider this line as an accu- 
rate and truly isotropic indicator of AGN activity, even 
if it could also originate from high excitation starburst 
emission or in shocks in low-metallicty starbursts (Lutz 
et al 1998). 

The [NeIII]15.55/zm line (Ne+ and Ne++ have ion- 
ization potentials of 14eV less than 0++ and 0+++ 
respectively) can be excited both from AGN activity 
(Gorjian, Cleary & Werner 2007) and from starbursts 
(Thornley et al. 2000). Superimposed on the AGN 
spectra, the lines of [NeII]12.81/xm, [SIII]18.71/xm and 
33.48/im and [SiII]34.82/im originate in gas with mod- 
erate ionization and most of their emission is generated 
by young newly formed stars, even if some contribution 
from the AGN is also expected (Spinoglio & Malkan 

1992) . The relationship between the [OIV]25.88/zm, 
[NeIII]15.55/im and [NeII]12.81^m lines in an eteroge- 
neous sample of Seyfert galaxies has been studied by 
Melendez et al. (2008a), who found that Seyfert l's 
and Seyfert 2's have different AGN and star formation 
contributions to the total emission. 

The interstellar medium produces the pure rotational 
lines of molecular hydrogen, as already shown by the 
early results of ISO spectroscopy in Gcnzel et al. (1998) 
and Rigopoulou et al. (2002), respectively. Wu et al. 
(2009) analyzed the IRS low resolution spectra of 103 
Seyfert galaxies from the 12MSG and measured the pol- 
icyclic aromatic hydrocarbons (hereafter PAH) emission 
features and the Silicate absorption strength. The PAHs 
have been proposed as star formation tracers by Puget 
& Leger (1989), while the Silicate absorption is sensitive 
to heavy dust obscuration of the nucleus. 

According to the simplest Unified Model for AGN, the 
Accreting Torus Model (ATM) (Malkan, Gorjian & Tarn 
1998), Seyfert 1 and 2 galaxies are the same kind of ob- 
jects, only viewed from different angles. The strongest 
demonstration of this is detection via optical spectropo- 
larimetry of the Broad Line Region emission-the defin- 
ing characteristic of Seyfert l'-in a significant minority of 
Seyfert 2 galaxies (Antonucci & Miller 1985; Antonucci 

1993) . A different scenario postulates an evolutionary 
difference: that Seyfert 2 are the early stages of the tran- 
sition of HII/Starburst galaxies into Seyfert l's. Two 
suggested evolutionary progressions are HII — ► Seyfert 2 

2 either currently or some time during the last several hundred 
years, because, even if the AGN ionizing continuum could be com- 
pletely switched off, the photoionized NLR could still be detected, 
due to its large extension and the long recombination time 



(Kauffmann et al. 2003; Storchi-Bergmann et al. 2001), 
or a fuller scenario of HII — > Seyfert 2 — ► Seyfert 1 (Hunt 
& Malkan 1999; Krongold, Dultzin-Hacyan & Marziani 
2002; Levenson, Weaver & Heckman 2001). Because the 
radiation due to the star formation processes is roughly 
isotropic, the ATM predicts no observational difference 
in the star formation tracers between Seyfert l's and 2's. 
If, on the other hand, star formation is stronger in Seyfert 
2's, as suggested in Buchanan et al. (2006), then some 
evolution from Seyfert 2's to Seyfert l's could be invoked. 
If general interstellar extinction towards the center of the 
galaxy is extremely high, optical data alone might not al- 
ways provide the correct (intrinsic) classification. 

Our sample is described in §2, the observations and the 
data reduction are briefly reported in §3, the direct re- 
sults of our observations, the estimates of the H 2 masses 
and temperatures and the measure of the continuum ex- 
tendedness are presented in §4, the diagnostic diagrams 
and the semianalytical models to interpret them are pre- 
sented and discussed in §5. In §6 the [NeV] is quanti- 
fied as an unambiguous AGN activity indicator, in §7 
we show that the mid-infrared diagnostics differentiate 
AGN l's from the other populations, and we derive the 
average spectra for each class of objects that can be used 
as templates also for predictions and comparisons with 
high redshift populations. Finally in §8 we present the 
line luminosity functions for our sample, and calculate 
the total accretion power generated in the local universe. 
The conclusions are summarized in §9. 

2. THE SEYFERT GALAXIES OF THE 12/iM GALAXY 
SAMPLE 

From the original Seyfert galaxies list of the RMS, 
we present 91 IRS high-resolution spectra, including one 
third of them which were published in Paper I. Our final 
sample is over 80% complete, large enough to give rea- 
sonable statistical results, with 41 Seyfert l's, 47 Seyfert 
2's and 3 galaxies which have been reclassified as optical 
starburst galaxies, according to NED 3 . 

Another improvement of this work is the classification 
that we adopt: we reclassify the Seyfert l's and 2's into 
more general "AGN l's" and "AGN 2's". We consider 
AGN l's to be all those with broad line regions, including 
those Seyfert 2's with hidden broad line regions (here- 
after HBLR), observed in polarized light. The remaining 
Seyfert 2's lacking any broad permitted lines, even in 
polarized light, are classified as AGN 2's. Our classifica- 
tion scheme is an attempt to identify a 'clean' category 
of intrinsically broad-line AGN. We follow Tran (2001) 
and Tran (2003), who made a spectropolarimctric sur- 
vey of the 12MSG. Out of the original 47 Seyfert 2's, they 
found HBLR in 19, 20 lacking an HBLR and they reclas- 
sify 11 objects as LINER, HII or starburst galaxies 4 . We 
classifiy as non-Sy for all these latter objects and similar 
ones (in the diagrams LINER and HII or starburst galax- 
ies will be distinguished). Exceptions are NGC1097 and 
MRK897, that we reclassify as being a AGN 1 and an 

3 NASA Extragalactic Database, IPAC, Caltech Pasadena, 
http:/ /nedwww. ipac.caltech.edu/ 

4 Among these If objects: three (Mrk897, NGC7496 and 
NGC7590) were presented in Paper I and the other eight are: 
NGCf056, NGC1097, NGC4922, NGC5005, NGC68fO, NGC7f30, 
MCG+0-29-23 and CGCG38f-05f, whose IRS spectra are pre- 
sented in this paper. 
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AGN 2 respectively. We refer to the Appendix A for the 
details on the classification of each one of these objects. 
Tran (2003) also adopt a somewhat arbitrary distinction 
between "bona fide" and "non-bona fide" Seyfert l's, 
based on the NED classification of Seyfert 1.8, Seyfert 
1.9 types and on radio loudness. For the remaining 13 
objects not considered as "bona fide" Seyfert 1 in Tran 
(2003), we prefer instead to carefully classify them on the 
basis of detection or not of optical broad lines, either in 
direct or polarized spectra. We present the details on the 
classification of these 13 galaxies in Appendix A, while we 
summarize here that we classify as AGN 1 seven Syl of 
the original RMS list (NGC526A, NGC1097, NGC1365, 
NGC2639, NGC7316, ES0545-G13 and ES0362-G18), 
we classify as a HBLR NGC5347 (Moran 2007) and 
we reclassify as AGN 2 one galaxy (NGC5506), while 
we consider non-BLR four objects for which there is no 
evidence of broad lines, but they lack polarization ob- 
servations (NGC1194, NGC4602, MRK1034 NED02 and 
MRK897) and finally we reclassify as non-Sy two galax- 
ies (NGC3511 and MRK1034 NED01). Although Tran 
(2003) distinguish the radio-loud 3C galaxies, we actu- 
ally classify 3C120, 3C234 and 3C445 as AGN 1, inde- 
pendently of their radio characteristics, because of the 
presence of broad line emission. 

We note that 7 Seyfert l's have Balmer lines with rel- 
atively small FWHM, under 2000 km/s but which are 
nonetheless produced in a BLR (Zhang & Wang 2006), 
These are usually classified as "narrow line Seyfert l's", 
but our sample contains too few to define another cate- 
gory. It will turn out that their infrared spectra do not 
appear different from those of normal Seyfert l's (see Sec. 
4). 

In summary our observed sample of " original" Seyfert 
galaxies contains 55 objects showing some evidence of 
broad lines, either in direct (34 Sy 1) or polarized light 
(21 HBLR) that we classify as AGN 1, 20 AGN 2 (non- 
HBLR) and 4 non-BLR not included in the AGN2 class. 
We will also consider in the following the 13 non-Sy 
galaxies, however these latter will not be used for any 
statistical derivation. 

Following the results of Wu et al. (2009), who have 
defined as "20 peakers" the Seyfert galaxies having 
a flux ratio F2o flm /F3Q^ m > 0.95, we have also identified 
these in our sample 5 to search for any difference with 
our classes of galaxies. We refer to Appendix B for the 
results on these objects. 

3. OBSERVATIONS AND DATA REDUCTION 

Most of the sample galaxies-the 29 galaxies presented 
in Paper I and the 23 galaxies in this paper- have been 
observed within the Spitzer Guaranteed Time Project 
30291 (P.I. Fazio). IRS high-resolution observations of 
another 37 objects, belonging to the 12MSG, were ex- 
tracted from the SSC (Spitzer Science Center) archive. 
For 25 of the latter objects, the observing mode was sim- 
ilar to the one of P30291, namely the off-source observa- 
tions have been collected to allow accurate background 
subtraction. For the remaining 12 sources from the 
archive (see Tab.l), no off-source observation was taken. 

5 The "20 (im peakers" in our sample arc: MRK335, MRK348, 
NGC424, NGC526A, MCG-02-08-039, F03450+0055, ESO033- 
G002, MRK0006, MRK704, MRK1239, 3C234, MCG-06-30-015, 
IC4329A, NGC6860, NGC7213 



For these objects, we give the line intensities, as for the 
other galaxies, but not the equivalent widths (EW) of 
the emission features and lines, nor do we present the 
photometric measurements at 19 and 19.5 /zm for com- 
puting the extendedness of the sources, because of their 
large uncertainty in the level of the continuum. MRK335 
was observed by the Program 50253 (P.I. Malkan), in the 
last cold cycle of Spitzer. 

The whole data reduction was done using the SMART 
packages 6 . We refer to Paper I for more details of the 
data reduction and analysis. 

4. OBSERVATIONAL RESULTS 

The journal of the observations of the 61 galaxies pre- 
sented here is shown in Tab.l, giving for each galaxy: the 
equatorial coordinates at the 2000 equinox; the redshift; 
the original RMS Seyfert class and the new classification; 
the IRAS fluxes at 12 and 25/zm; the observing date and 
the number of cycles and integration times per cycle. For 
the sources already presented in Paper I, we refer to the 
table 11 in the appendix for the adopted classification. 

The spectra of the new 61 galaxies are shown in Fig.l. 
For all galaxies for which the off-source spectra have 
been subtracted, both the SH and the LH spectra are 
presented. For the galaxies for which no off-source ob- 
servation was available, we show in the figure the on- 
source SH spectrum only. This latter is only marginally 
affected by the lack of background subtraction, because 
the theoretical background, as measured using the back- 
ground estimator provided by the SSC, is less than 10% 
of the total measured SH emission. For these galaxies, 
we show only the detected lines in the LH range; we do 
not present the whole spectrum, because it is affected 
by a higher level of background (estimated to be about 
20-30% of the total measured emission). The spectra of 
MRK335, F05563-G018, NGC5135, IC4329A, NGC5347 
and NGC5506 show only the LH detected lines, because 
after background subtraction, the LH orders are not well 
inter-calibrated and thus the continuum cannot be de- 
fined properly. 

Tab. 2 reports the fluxes of the fine structure lines, mea- 
sured with a Gaussian fit, for both the SH and LH spec- 
tra. Tab. 3 gives the fluxes of the molecular H2 rotational 
lines S(0), S(l), S(2) and S(3) and the PAH 11.25/zm in- 
tegrated flux, measured with a moment fit, and its equiv- 
alent width. We consider as detections the measurements 
with a signal to noise ratio higher than 3. We also report 
in Tab. 11 in the appendix the PAH 11.25/im integrated 
flux and equivalent width as well as the H2 masses de- 
rived from the H2 rotational lines observations for the 30 
spectra of the galaxies presented in Paper I. We have re- 
peated the measurements of fluxes and equivalent widths 
of the PAH 11.25/im to be consistent with the following 
method adopted here, and also because we found that 
some of the equivalent widths given in Paper I were af- 
fected by non-traceable errors. 

We measured the PAH fluxes by removing from the 
spectra the continuum under a baseline traced from the 
continuum shortwards of the PAH feature to the contin- 
uum longwards of the [Nell] line. Such a large interval 

6 SMART is available on the SSC website and developed by the 
Infrared Spectrograph (IRS) Instrument Team at Cornell Univer- 
sity (Higdon et al. 2004) 
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has been chosen because, in addition to the feature at 
11.25pm, two other PAH features are present (approxi- 
mately at 12.0pm and at 12.5pm) and they increase the 
level of the apparent continuum under the 11.25pm fea- 
ture. The integration range of the PAH emission depends 
on the feature's intensity, for the brightest sources it can 
be as wide as 0.5 pm (11.15pm - 11.65pm). Chosing this 
large baseline allowed us to avoid the other PAH con- 
tributions and remove the correct continuum from the 
galaxy. By dividing the resulting PAH integrated flux by 
the continuum flux density at the midpoint wavelength 
of the baseline, we can obtain the EW of the feature. 

4.1. H 2 excitation diagrams: H 2 temperatures and 
masses 

Using the H2 rotational lines intensities, we can esti- 
mate the temperature and the mass of the H2 line emit- 
ting regions (see Paper I for the details of the method 
used). The derivation of the masses presented in Paper 
I was affected by some errors, therefore we recomputed 
their values, and list them in Table 9. 

In Fig. 2 the excitation diagrams of each of the sources 
are shown. Tab. 4 presents the derived temperatures T(3- 
2), T(2-l), T(l-O) of the gas in those regions where the 
rotational lines S(3), S(2) and S(l) are emitted, respec- 
tively. The 0-0 S(l) is the most intense H2 line in every 
source, and so we estimated the mass of the molecular 
regions from its intensity. We calculated the masses from 
the Boltzmann equation using the two transitions of S(0) 
and S(l), listing the results in Tab. 4. Therefore we ob- 
tain the masses only for the 50 galaxies where both these 
rotational lines were detected, while we give an upper 
limit to the mass when only the S(l) line was detected. 
The temperature values that we derive here are compa- 
rable with the ones we report in Paper I. 

The estimated H 2 masses are of the same order of mag- 
nitude as those calculated by Rigopoulou ct al. (2002) 
in Seyfert and starburst galaxies and by Higdon et al. 
(2006) in ULIRGs, ranging from 10 7 to 10 9 M . The 
average mass (in units of 10 8 M ) for the AGN 1 class is 
1.8±1.3, for AGN 2 is 0.65±0.47, excluding from the av- 
erage the outlier NGC1142 because it has a mass an order 
of magnitude larger than the others. For non-Sy galaxies 
is not possible compile a sensible average, because their 
masses spread over a wide range of values; MRK1034 
NED1 has a mass of 5.36, CGCG381-051 3.8, NCG3511 
0.15 and NCG7590 0.12. For the two Seyfert galaxies in 
common with Rigopoulou et al. (2002), NGC1365 and 
NGC7582, our results are in complete agreement, even 
though they used ISO spectra, for which the H2 S(l) data 
were obtained throught an aperture of 14"x27", that is 
7 times larger than the aperture of SH. This implies that 
the H2 emission is concentrated in the inner ~ 50" 2 of 
these two large (several arcminutes in diameter) galaxies. 

In figure 3a, b we present the H 2 flux vs PAH flux and 
the H2 luminosity vs PAH luminosity. We confirm the 
results of Paper I with our new classification: there are no 
differing trends which could discriminate between AGN 
1 and AGN 2. 

4.2. Source extendedness 

The SH and LH spectra overlap in the range 17-19pm, 
allowing us to form an estimate of the extendedness of 
the sources called R, the ratio of the flux measured in 



LH to the flux measured in SH in an adjacent spectral 
interval. This parameter can be defined only for sources 
after an appropriate background subtraction. We refer 
to Paper I for the details. The flux densities measured 
at 19pm in SH and those at 19.5pm in LH are reported 
in Tab. 5, with the R ratio and the corresponding ex- 
tendedness class (see Paper I). We report the average 
extendedness factors in the classification table (Tab. 6), 
for each class of galaxies and also for the group of " 20pm 
peakers". Both SH and LH spectra were already cor- 
rected for light lost outside the slits, assuming a point 
source: this correction takes into account the increasing 
size of the PSF as wavelength increases from the blue to 
the red end of each spectrum. In order to better under- 
stand the validity of the spectroscopic extendedness ra- 
tio R, we analyzed Spitzer direct images of four compact 
(R~l) sources taken with IRAC at 8pm and MIPS at 
24pm. We measured in these images the fluxes through 
the same apertures as observed by IRS in SH and in LH, 
flux-corrected with the prescriptions of the SSC for point- 
sources. The mean extendedness of these four sources 
measured through IRS is < R >irs= 1-00 ± 0.01, while 
for IRAC at 8pm and MIPS at 24pm are respectively 
< R >irac= 1.00 ±0.03 and < R >mips= 1.02 ±0.06. 
Therefore the IRS estimate of the extendedness is reliable 
for compact sources. We have also attempted an analo- 
gous test on the extended sources, although this estimate 
has serious limitations because the aperture corrections 
to photometry of extended sources are not well defined 
for IRAC 7 and MIPS. We find that, using both IRAC 
at 8pm and MIPS at 24pm, the extendedness classes we 
have defined are reproduced for our sources, although 
there are some differences in individual values. This im- 
plies a broad agreement among the IRS, IRAC at 8pm 
and MIPS at 24pm measurements of the extendedness. 
A paper with the analysis of the IRAC four channels 
and MIPS at 24pm images of the Seyfert galaxies of our 
sample is in preparation. 

5. AGN DIAGNOSTIC DIAGRAMS: DATA AND MODELS 

One of our aims is to develop a method to disentan- 
gle the contributions of the AGN and the Starburst to 
the total infrared emission of the Seyfert galaxies of our 
sample using mid-infrared spectral features. In this sec- 
tion we use the diagnostic diagrams, together with semi- 
analytic models, to estimate for the sample galaxies the 
AGN contribution, e.g. through the [NcV] to [Nell] line 
ratio, and the star formation contribution, e.g. with the 
equivalent widths of the 11.25pm PAH feature and of the 
[Nell] emission line, and to search for differences between 
AGN l's compared with AGN 2's. 

To estimate the percentage of the contribution of 
the AGN and the Starburst to the observed emis- 
sion of a galaxy at 19pm, we computed a mathe- 
matically simple analytic model to describe both the 
emission due to the AGN and that due the star for- 
mation in each of the following observed quantities: 
the extendedness of the source, the equivalent widths 
of the PAH at 11.25pm and of the [Nell] line at 
12.81pm, the line ratios [NeV]14.32pm/[NeII]12.81pm 
and [OIV]25.89pm/[NeII]12.81pm and the spectral in- 
dex a at (60-25)pm. We constructed analytic models for 

7 For IRAC, sec: http://ssc.spitzer.caltech.edu/irac/calib/cxtcal/index.html 



Spitzer spectra of Seyfert galaxies 



5 



these quantities, because they provide the best separa- 
tion of emission from AGN and non-AGN, and therefore 
the best estimates of the AGN percentage contribution 
to the total mid-infrared emission at 19/xm. The simple 
equations for each of the models are given in the Ap- 
pendix C. In the following plots of those quantities, we 
compare the observations with the semi-analytic models, 
which are plotted with solid lines. 

5.1. Line ratios versus PAH equivalent widths 

In Fig.4a the line ratio [NcV]14.32/xm/[NeII]12.82^m 
is shown MS iX function of the equivalent width of the 
PAH feature at 11.25/xm. This ratio is the best AGN 
tracer in the IRS wavelength range, because [NeV] can 
be excited only by the AGN ionizing continuum. Its ratio 
to [NeII]12.82/im is not directly affected by abundances. 
76% of the AGN l's have an absolute value of the EW of 
the PAH at 11.25/xm of \EWPAH\ < 0.1 Vn and a neon 
ratio [NeV] /[Nell] > 0.15. On the other hand, 87% of 
non-Sy and non-bona fide Syl have \EWPAH \ > 0.14/um 
and [NeV] /[Nell] < 0.5. AGN 2's (i.e. non-HBLR) show 
a wide range of both the EW of the PAH and of the neon 
line ratio. Fig. 4b shows a similar diagram with the PAH 
EW vs [OIV]/[NeII] line ratio. This diagram presents the 
same characteristics of the former, but with fewer upper 
limits, confirming its results. There is a Syl, NGC7213, 
that has not been detected in the [OIV] line and has been 
reported to be at an intermediate stage between LINER 
and Seyfert (Starling et al. 2005). 

Both the diagrams of Fig. 4 show semi-analytic mod- 
els which reproduce the empirical data and estimate the 
percentage of the AGN and the starburst contributions 
to the total infrared emission at 19[im. If we define 
as an Seyfert- dominated galaxy any of the sources with 
an AGN contribution equal or more than 50% (see the 
dashed lines in the diagrams), this corresponds to a the 
line ratios [NeV] /[Nell] > 0.054 and [OIV] /[Nell] > 0.28. 
Using this simple argument, we confirm that all the AGN 
2's (non-HBLR's) arc Seyfert-dominated, even if they 
span a wide range of AGN percentage in the diagnos- 
tic diagrams. We report in the Table 6 the average PAH 
EW and [NeV] /[Nell] ratio for each class of galaxies. 

Also the diagrams PAH EW vs 
[NeV]14.32/zrn/[SiII]34.8/zm and PAH EW vs 
[NeIII]15.55/xm/[NeII]12.82/xm, presented in Fig.5a,b, 
show similar characteristics as the previous two: the line 
ratio increases when the PAH EW decreases. All these 
diagrams (Figs. 4, 5) show the general inverse relation 
between the AGN dominance, as measured from the 
ionization sensitive line ratios, and the star formation 
dominance, estimated from the PAH EW. 

5.2. Line ratios as density indicators 

In the IRS high-resolution spectral range, there are 
two fine structure line doublets, from which the electron 
density can be derived. These are [NeV] at 14.32/im and 
24.31^m and [SIII] at 18.71/xm and 33.48^m. Follow- 
ing Dudik et al. (2007), as already reported in paper 
I, we present the diagram of the neon doublet ratio ver- 
sus the sulfur doublet ratio in Fig. 6. In the diagram are 
shown the low-density limits for both line ratios and the 
estimated electron density as calculated by Dudik et al. 
(2007) assuming a temperature T = 10 4 K. The sulfur 
doublet ratio has been computed correcting the flux of 



the [SIII]18.71/im taken throught the smaller aperture 
SH, multiplying it by its extendedness factor (see Ta- 
ble 5) to be compared with the larger aperture LH. If 
the ratio of two lines of the same species is below the 
low density limit, either there is extinction that prefer- 
entially reduces the flux of the shorter wavelength tran- 
sition, and/or the assumed temperature is wrong, and 
these doublet ratios could not be used to measure the 
electron density. 

We find that 34% of the AGNl's lie under the neon ra- 
tio low density limit and 15% under the sulfur low density 
limit, and that 30% of the AGN2's lie under the sulfur 
limit and 30% under the neon's. Our current sample is 
much larger than that of Paper 1, making these results 
statistically significant. 

Taken at face value, these diagrams indicate low elec- 
tron densities (lO*- 3 ^ cm -3 ) in the [NeV] emitting re- 
gion. For a typical ionizing luminosity of 10 43 erg s _1 and 
an ionization parameter of 10~ 2 , this implies a character- 
istic radius of about 100 parsecs. Since the recombina- 
tion time is long, for densities of lO*- 3-4 ' cm~ 3 it ranges 
from 10-100 years 8 , 'fossil' NLR emission will continue 
to be detectable in AGN 2 which could have 'shut down' 
production of ionizing photons for the last 300 years. We 
suspect this may be the explanation of many of the 'pure 
Seyfert 2' members of our AGN 2 class. 

Thus a significant percentage of each of the Seyfert 
types falls under the low-density limit for either one or 
the other doublet, therefore, at least for the galaxies 
whose line ratio(s) are below the limit, we cannot es- 
timate correctly the electron densities. However all the 
objects are below the low density limit only in one line 
ratio and not in the other line ratio. If extinction were re- 
sponsible, both ratios would show this effect. We do not 
have any explanation for this unexpected behaviour and 
the possibility that some currently adopted atomic pa- 
rameters for these lines (transition probabilities and/or 
collision strengths) might be inaccurate must be consid- 
ered. 

In Fig. 7a we show the [SIII] doublet ratio as a func- 
tion of the [SIV]10.51/im/[SiII]34.8^m ratio, with the 
[SIII]18.71/im flux corrected for the extendedness factor 
(cf Sec. 4.2). These lines are produced either in the 
narrow-line regions of the active nuclei or in the HII re- 
gions. In Fig. 7b the [NeV] doublet ratio as function of 
the [NeV]14.32^m/[NeII]12.82^m line ratio is given. Fig- 
ures 7a and b show no differences in the density-sensitive 
line ratios between the different Seyfert populations. 

5.3. Ionization Diagrams 

[NcV]14.32^m and [OIV]25.89^m are the best AGN 
tracers in the high-resolution IRS spectra. We normal- 
ized them both to the [NeII]12.81/xm, that is mainly pro- 
duced in the HH/star- forming regions, to produce two 
ionization-sensitive ratios (in Fig. 8a). All the sources lie 
along the same sequence. As expected, the AGN l's oc- 
cupy the highest ionization region of the diagram as in 
the other diagrams. Most AGN 2's lie in the same region 
of the AGN l's. This diagram can be used to estimate 
the ionizing power in the NLR, both in AGN 1 and AGN 
2 objects. As mentioned above, we computed a semi- 

A typical recombination time is r r ~ 10 5 /n e years (Osterbrock 
& Ferland 2006) 
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analytic model to reproduce the data in this diagram. 
We choose as the threshold to consider a source as an 
Seyfert-dominated AGN again a 50% AGN contribution 
to the 19/jm continuum, and this includes all the AGN l's 
and all the detected AGN 2's. A clear implication is that 
using optical spectroscopy to identify Seyfcrt nuclei, as 
we have done in the 12MGS, will systematically miss out 
a significant population of non-HBLR Seyfcrt 2's (weak 
AGN 2's) which have less than 50% AGN contributions 
to their infrared emission at 19/jm. 

Vcilleux et al. (2009) present the same diagram for 
their sample of ULIRGs belonging to the QUEST survey, 
and compute a model to estimate the AGN and starburst 
contributions to the total emission. Their values of the 
line ratios are in the same range as ours. However, the 
percentage of AGN contribution they estimate is on av- 
erage lower than the percentage we compute. This differ- 
ence is likely due to the different populations of the sam- 
ples. In fact ULIRGs have a stronger starburst compo- 
nent than Seyfert's and/or more heavily obscured AGN 
emission, even in the mid-infrared. 

5.4. Line Equivalent Width Diagnostics 

The line equivalent widths can also be used to esti- 
mate the contributions of the AGN and the starburst 
to the total galaxy emission. Since the original IRAS 
studies of Seyfert galaxies (Spinoglio & Malkan 1989; 
Rush, Malkan & Spinoglio 1993; Spinoglio et al. 1995; 
Spinoglio, Andreani & Malkan 2002), we have known 
that the mid-IR continuum can be dominated by the 
AGN continuum (e.g. from 12/jm to 25/xm), while at 
longer wavelengths the continuum due to the galactic 
disk and reprocessed thermal emission from dust around 
the young stellar populations increases, and it dominates 
the total emission at wavelenghts of 60-100/jm. There- 
fore the equivalent widths of emission lines in the mid- 
IR are differently affected by the underlying continuum, 
as a function of the wavelength: around 12/im the EW 
will be more depressed by the strong AGN continuum, 
than they would be around 25/^m. The lines originat- 
ing in the galactic disk and its star- forming regions, such 
as [NeII]12.8/zm and the H 2 S(l) at 17.02/im, will have 
a smaller equivalent width in objects with a stronger 
AGN continuum, while the lines tracing the AGN, like 
[NeV]14.32/jm and the [OIV]25.89/jm will not suffer this 
depression because they are enhanced in proportion to 
the strength itself of the AGN, keeping their equivalent 
widths more constant. 

In our analysis of the equivalent widths, we con- 
sider only objects for which an accurate measurement 
of the continuumm was possible through background 
subtraction and with detected lines, neglecting the 
ones with upper limits. The equivalent width of the 
[NeII]12.8lMTO, the [NeV]14.32/jm, the [NeIII]15.55/zm, 
the [OIV]25.89/jto and H2 17.02/jto are reported in the 
Tab. 7. 

We present in Fig. 8b the PAH EW as a function of the 
[NeII]12.82/im equivalent width. This diagram is able 
to disentangle Seyfert-dominated AGN from starbursts. 
Comparing it with our semi-analytic model, we find that 
almost all the sources lie in the region of AGN contri- 
bution greater than 50%. We further find all the AGN 
l's and 50% of the AGN 2's lie in the region with the 
modelled AGN contribution >75%, i.e. for |EW [NeII]| 



< 0.08. The likely explanation is that both the HII 
regions/star- formation tracers- [Nell] and PAH-happen 
to have wavelengths in a continuum range that can be 
strongly enhanced by the presence of an AGN contri- 
bution. The brighter the nuclear AGN continuum, the 
smaller the equivalent width of the spectral feature, as 
can be seen in the average spectra of the different popula- 
tions of Fig. 14 in §7.1. We report in table 6 the [NeII]EW 
for each class of galaxies. 

We computed a mean PAH EW at 11.25/jto for the 
AGN 1 of < |PAH EW| > = 0.17±0.20 (Seyfert 1: 
0.17±0.21 and HBLR: 0.16±0.19), and for the AGN 2 

< |PAH EW| > = 0.34±0.28. Wu et al. (2009) find, 
for the 55 objects that are in common, a similar mean 
PAH EW at 11.25/jto. From their data, we calculated for 
the AGN 1 class < |PAH EW| > = 0.19±0.20 (Seyfert 
1: 0.20±0.20 and HBLR: 0.17±0.21), and for the AGN 2 

< |PAH EW| > = 0.43±0.27. For the detailed compar- 
ison between ours and their data, we refer to Appendix 
D. 

Using a Kolmogorov-Smirnov test, we derive that the 
probability that the two kind of AGN l's- the Seyfert l's 
and the HBLR-Seyfert 2's- belong to the same family is 
78%. Applying the same test to the Sy 1 and the AGN 
2 or to the HBLR-Sy2's and the AGN 2, results in both 
cases in a probability of less then 1% that the two groups 
are drawn from the same population. This reinforces our 
hypothesis that the HBLR belong with the Syl, not the 
other Sy2s. We will further discuss this topic in Par. 7. 

5.4.1. Line equivalent width, extendedness factor and 
far-IR slope 

We present in Figures 9-10 the extendedness versus 
the equivalent widths of 11.25/im PAH, [NeII]12.8/jm, 
and [NeIII]15.5/zm, respectively. The galaxies with the 
largest extendedness, which we define as class III, namely 
NGC1056, NGC1142, NGC1667 and NGC5005, are the 
ones with the largest equivalent widths of the PAH at 
11.25/xm and of the [NeII]12.8/im. These emission fea- 
tures are generated principally in the interstellar medium 
and in the HII regions of the galaxies, but not much from 
their active nuclei. Among these four galaxies, NGC1056 
and NGC5005 have no detected [NeV] lines, and are 
reclassified as non-Sy (see tables 1 and 2); NGC1142, 
NGC1667 are AGN 2. Including the galaxies of Paper I, 
among the galaxies of class III there is also a non-BLR, 
NGC4602, and a Syl, NGC4748. We find the same trend 
as in Paper I, that with increasing size of the 19^m emit- 
ting region, both the PAH and [Nell] equivalent widths 
increase. That is, in these diagrams AGN l's occupy the 
region where the point sources lie. In contrast the non- 
Sy 's with a measured extendedness, occupy the region of 
significantly extended sources with large PAH or [Nell] 
equivalent widths. This separation works properly in the 
diagrams of extendedness vs PAH EW and extendedness 
vs [Nell] EW, shown in Figg.9a and b. In fact just 1 of 23 
Sy l's and 2 out of 17 HBLR's do not lie in the "AGN 1 
region", with extendedness <1.3. But 2/3 of the non-Sy 
have extendedness >1.3. AGN 2 split almost equally in 
the two regions: 8 of 18 lie in the AGN 1 region and 10 
in the other. 

Comparing the diagrams of the extendedness versus 
EW PAH and versus EW [Nell] with our semi-analytic 
model, we find that these diagrams do separate the 
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Seyfert-dominated AGN from Starburst-dominated ones: 
a modelled AGN contribution greater than 75% corre- 
sponds to |EW PAH| <0.4//m and |EW [NeII]| <0.08^m. 

The diagram of [Nelll] equivalent widths versus the 
19/zm source extendedness (Fig. 9) shows a similar trend 
of the figure with [Nell], but with a much higher scatter. 
On the contrary, the diagrams of the [NeV] and [OIV] 
equivalent widths versus the 19/im source extendedness 
(that we do not show) do not show any clear trend. This 
confirms that the high-ionization lines are almost totally 
produced in the active nuclei with no significant contribu- 
tion from in extended regions. [Nelll] is an intermediate 
case, with contributions from both the Seyfert nucleus 
and also the extended star- forming regions. 

The IRAS 12/im luminosity and the 
[NeV]14.32/zm/[NeII]12.82^m line ratio are not corre- 
lated for any Seyfert type (see Fig. 11), confirming the 
finding of Paper I. We notice that high-luminosity 
objects, e.g. those with £i2 M m > lO^ergs^ 1 , 
have the ratio [NeV] 14.32/xm/ [Nell] 12.82^m > 
0.2, i.e. are AGN-dominated. Fig. 12a, b shows 
the diagrams of [NeV]14.32^m/[NeII]12.82/zm and 
[NeIII]15.55/zm/[NeII]12.82/im vs (60-25)/xm spectral 
index. In both diagrams AGN l's lie in the upper right 
region, with [NeV]/[NeII] > 0.05 and [NeIII]/[NeII] > 
0.1, and the spectral index ck6o-25 > —2.2. The non-Sy's 
cluster in a region [NeV] /[Nell] < 0.2 and [NeIII]/[NeII] 
< 0.5, with a steeper spectral index, < —1.2. 

AGN2 lie across both the regions. The average spectral 
indices are given in table 6. 

5.5. The AGN and starburst contribution to the 19\im, 

flux 

By inverting the semi-analytic models, i.e. solving 
the analytical expressions with the true values of the 
observed quantities, we obtained a value of 5ft, defined 
as starburst to AGN continuum ratio at 19/xm (see Ap- 
pendix C) for each of the five observed quantities of ev- 
ery source. We then computed the mean value of 5ft to 
estimate the relative percentage of AGN and Starburst 
emissions. From the value of 5ft we computed in Tab. 8 the 
percentages of the Starburst and of the AGN emission. 
When, for a particular source, the scatter is greater than 
half of the mean and is due to a single 5ft value of one of 
the models, which does not describe that source because 
of a noisy detection, we removed that discrepant value 
of 5ft and recomputed the mean. The sample of sources 
to which we apply this analysis is reduced from the ob- 
served sample of 89 because the models depend on the 
extendedness and the PAH equivalent width, which are 
not always detected (for the PAH) or measurable (for the 
extendedness) . We found that the model can disentangle 
the AGN and the Starburst emission for 31 AGN1 (17 
bona fide Syl's and 14 HBLR), 3 non-BLR, 15 AGN 2's 
and 9 objects re-classified as non-Seyferts. As can be seen 
from the histogram in Fig. 13 Syl's have a mean AGN 
contribution at 19^to of 92% ± 6%; HBLR's 92% ± 8%; 
AGN's 79% ± 16%, non-Sy's 69% ± 16%, and non-BLR's 
62%±7%. These average percentages are also given in ta- 
ble 6. It is perhaps an uncomfortable surprise that up to 
halfoi the 19^m continuum in the IRS slit can come from 
a normal AGN, but still not be strong enough to make it 
unambiguously classified as an AGN from spectroscopy. 
As discussed below, this has lead to substantial confusion 



in the literature when Seyfert "AGN" are discussed and 
compared, using different selection observations. 

The differences that we find between the AGN 1 and 
the AGN 2, (even if the latter cannot be considered from 
our data as an homogeneous population), i.e. the lower 
ionization ratios, the increased PAH and [Nell] EQW, 
the extendedness of the 19/im emission, can all be re- 
lated to a weaker strength of the AGN component with 
respect to the starburst component. As a matter of fact, 
all our AGN 2's have evidence of an AGN both in the 
optical (as seen in their optical line ratios in the BPT di- 
agrams (Baldwin, Phillips & Terlevich 1981), where we 
have used the optical spectral observations of Rodriguez, 
Malkan & Spinoglio (2009)) and also in the hard X- 
rays, at 2-10keV (Shu et al. 2007; Malizia et al. 2007; 
Bianchi et al. 2005), 2-8keV (Cardamone, Moran & Kay 
2007), 15-136keV (Deluit 2004) or 40-100keV (Bird et 
al. 2008). 

5.6. Line equivalent widths for the different AGN types 

The average equivalent width of the fine structure lines 
in the various AGN classes will change depending on 
whether the line emission, and also the underlying con- 
tinuum, is dominated by the nonstellar nucleus or the 
hot stars. Thus the equivalent width of [NeV]14.32/im 
is the same for AGN l's: 0.020±0.021 and AGN 2's: 
0.025±0.010/im, because both the line and the under- 
lying continuum emission are proportional to the lumi- 
nosity of the AGN component. The EW drops in non- 
Seyfert galaxies because of the complete cutoff of [NeV] 
emission-it is detected in two non-Seyfert galaxies which 
are classified as LINER (see Appendix A), but measured 
only in NGC7130, for which we have an off-source mea- 
surement. 

The [OIV] EW decreases somewhat because the un- 
derlying continuum has a larger starburst contribu- 
tion passing from AGN1 to non-Sy, in fact the [OIV] 
EW of the AGN1 is 0.116±0.134^m, of the AGN2 
is 0.093±0.054^m and of the non-Seyfert galaxies: 
0.033±0.035^rn. In contrast, the [Nell] and the H 2 
equivalent widths increase for the same sequence of ob- 
jects, that is [Nell] EW for the AGN1 is 0.028±0.026^m, 
for the AGN2 it is 0.081±0.057^m, for the non-Seyfert 
galaxies it is 0.116±0.029^m and the average EQW 
H 2 is a factor ten greater in the non-HBLR's and the 
normal galaxies than in Syl's and HBLR's (AGN1: 
0.009±0.008^to; AGN2: 0.046±0.069/xm; normal galax- 
ies: 0.058±0.062^to). These trends are summarized in 
the 'classification table' given in Table 6. 

6. [NEV] AS INDICATOR OF AGN ACTIVITY 

Because of the very high ionization potential of NeV, 
we consider it as the best emission feature to distinguish 
active galaxies from starburst galaxies. Strong starburst 
activity could possibly excite some high-ionization lines 
such as [OIV]25.89^m (Lutz et al 1998), but the ion- 
izing spectrum of O and B stars is not hard enough to 
produce much Ne 4+ . Shocks would need to have excep- 
tionally high velocities. Among the 91 sources we ana- 
lyzed here and in Paper I classified as Seyfert galaxies, 16 
objects have no [NeV] emission. As already mentioned 
in Sec.2, 10 (NGC1056, MCG+00-23-029, NGC5005, 
NGC6810, CGCG381-051, MCG-03-34-063, NCG7496, 
NCG7590,NGC3511 and MRK1034 NED01) of these 16 
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galaxies have already been reclassified as LINER, HII or 
starburst galaxies by Tran (2001) and Tran (2003). 4 
are Syl's: NGC1097 and NGC2639 lie in the HII region 
area of the BPT diagram based on [NII]6584A/H a vs 
[OIII]5007A/H/3 (Baldwin, Phillips & Terlevich 1981); 
NGC7213, which also lacks the [OIV] line detection and 
has been classified as intermediate between LINER and 
Seyfert (Starling et al. 2005). The other undetected 
sources are the AGN 2 NGC4501 and MRK938, which 
also lack a detectable [OIV] line. 

In contrast, two of the non-Sy's show [NeV]14.32/im 
emission: NGC4922 and NGC7130 are LINERs, whose 
active nuclei can produce highly photo-ionized gas emis- 
sion (Dudik, Satyapal & Marcu 2009). There are also 
4 sources with detections of [NeV] at 14.32/im but not 
24.31/im. Nevertheless we consider them as [NeV] emit- 
ters, even if not detected at 24.31/im, because the noise 
of LH is greater than the noise of SH. 

In conclusion, we can consider the [NeV] emission lines 
as a strong indicator for a galaxy to be classified as an 
AGN. In fact it is detected in 88% of the AGN l's, 90% of 
the AGN 2's, only the 17% of the non-Sy's. This means 
that deep spectroscopic searches, e.g., for [NeV]14^m, 
can discover relatively weaker AGN with lower luminosi- 
ties, which only produce less than 45% of the total 19^m 
emission observed. 

Thus our infrared classification of galactic nuclei turns 
out to be in very close agreement to the classifications 
of them originally made from optical spectroscopy. The 
12MGS AGN sample, which is originally defined based 
on optical spectra, docs not contain very heavily ob- 
scured ("buried") Seyfert nuclei, more or less by defi- 
nition. Deep searches for [NeV]14/xm emission in the 
remaining "non-Seyfert" members of 12MGS would be 
required to find out how common buried Seyferts are. A 
first step has recently been taken by (Goulding & Alexan- 
der 2009), who uncovered a significant fraction of [NeV]- 
emitting luminous infrared galaxies whose Seyfert nuclei 
had been missed by previous optical spectroscopy, due to 
lack of sensitivity and heavy dust reddening. 

7. AGN 1 STATISTICS 

We give in Tab. 9 the probabilities P that the two 
AGN1 sub-populations, the Syl's and the HBLR Sy2's, 
do not differ significantly from one another, for each of 
the observed quantities. To derive those probabilicties 
we used a Kolmogorov-Smirnov test, which calculates 
the probability of two sets of data values arising from 
the same intrinsic distribution. The higher the probabil- 
ity, the closer are the two sets of data. We have grouped 
similar diagrams together to derive an average proba- 
bility < P >, with its standard deviation. For the first 
group of the relations combining fine structure line ratios 
with PAH equivalent width (EW PAH vs [NeV] / [Nell] , 
vs [NeIII]/[NeII], vs [OIV] / [Nell] , vs [NeV] /[Sill]), we 
obtain < P > = 73%±31% that the AGN1 populations- 
Sy l's and HBLR Sy2's-belong in the same class. For the 
second group, relations of the extendedness vs the equiv- 
alent width of the [Nell], [Nelll], [OIV], PAH we find 

< P > = 71%±12%. For the the relations between the 
H 2 S(l) line and the PAH flux and luminosity we obtain 

< P > = 77%±2%. Grouping together the relations be- 
tween the neon ratios ([NeV] /[Nell] and [NeIII]/[NeII]) 
vs the spectral index at (60-25)/im, gives < P > = 



99%±2%. A probability of P = 78% is obtained using 
the diagrams EW [Nell] vs EW PAH and ( [NeV] / [Nell] . 
The high probabilities obtained strengthen our classifica- 
tion which joins both Sy l's and HBLR's as members of 
a single AGN 1 category. We computed the same statis- 
tical tests to calculate the probability that the Syl and 
AGN2 or HBLR and AGN2 objects belong to the same 
family. Both of these probabilities are on average lower. 
The indistinguishability between Sy l's and HBLR has 
been also verified by plotting the average of the spec- 
tra of all AGNl's with an AGN percentage greater than 
95%, and the average spectrum of those AGNl's with an 
AGN percentage below that value. The resulting plot 
in Figure 14 looks almost identical to the plot with the 
average spectra of Sy l's and HBLR's. 

7.1. Average Spectra for each class 

To study representative spectral features of all classes 
of galaxies in our sample, and provide templates useful 
for comparisons and predictions of high-redshift galaxy 
populations, we computed the average spectrum of each 
of the classes normalized to the flux at 27/im. All the nor- 
malized average spectra are presented together in Fig. 14, 
and are available electronically. They are compared to 
the average spectrum of a sample of starburst galaxies 
from the IRS high resolution data by Bernard-Salas et 
al. (2009). 

The slopes of the continua of the average spectra are 
steepest in starburst and non Seyfert's and in AGN 2, 
with spectral indices 9 (10-35)^m of -2.88, -2.25 and - 
1.65, respectively. We cannot fit the continuum of the 
bona fide Syl's and the HBLR with a single power-law. 
Instead we must distinguish the (10-18)/^m and the (18- 
35)/xm continuum slopes. The shorter wavelength range 
in the bona fide Syl's has a spectral index of the averaged 
spectra of -1.21 and for the HBLR of -1.95, while the 
longer wavelength range has continuum spectral indexes 
of -0.65 and -0.90, respectively. 

Wu et al. (2009) find an average spectral index of the 
Seyfert l's to be -0.85±0.61 and of the Seyfert 2's to be 
-1.53±0.84, quite similar to the original findings based 
on IRAS data of Edelson & Malkan (1986). 

Fig. 14 shows a clear sequence. The highly ionized lines, 
such as [SIV], [NeV] and [OIV] are intense in the mean 
spectra of the Seyfert l's, HBLR Sy 2's, and AGN 2s. 
The PAH feature is stronger in starburst galaxies, non 
Seyfert's and AGN 2, and is weaker in AGN 1. Thus, go- 
ing from Seyfert-dominated to starburst-dominated ob- 
jects, the continua steepen from AGN 1 to AGN 2 to 
non-Seyfert's and starbursts, the higher ionization lines 
decrease not only in their flux, but also in their equiv- 
alent widths (see §5.6), while the PAH feature remains 
almost constant in flux, while its equivalent width in- 
creases. It seems therefore that the PAH grains do sur- 
vive in the highly ionized medium of AGNs, however they 
appear weaker in the most powerful ones, because they 
are masked by the strong underlying AGN continuum. 
The H2 S(l) line has a behavior similar to the PAH, 
and in fact its equivalent width increases with decreas- 
ing AGN activity (see §5.6) . The presence of a BLR is 
closely correlated with strong thermal continuum from 
hot dust grains which emit around 10/xm. 

9 We define as spectral index a = ■ 
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7.2. Line Luminosity Functions 

From our statistically complete sample of galaxies, we 
are able to compute the first infrared line luminosity 
functions for Seyfert l's and Seyfert 2's. The 12MGS 
includes 53 Seyfert l's and 63 Seyfert 2's, but only 42 
Seyfert l's and 50 Seyfert 2's have been observed by 
IRS, therefore we have to correct the luminosity function 
for completeness. Tab. 10 gives the luminosity functions 
for Seyfert l's and Seyfert 2's, for several of the bright 
emission lines. In Figg. 15-20 the luminosity functions of 
the [SIV], [Nell], [NeV]14.32^m, [Nelll], [SIII]18.71/zm, 
[NeV]24.31^m, [OIV], [SIII]33.5/mi, [Sill], H 2 17.04/^m 
and PAH 11.25/zm are shown. All the fine structure lines 
have a luminosity in the range 10 39 -10 43 erg sec -1 and 
the corresponding space densities range from 2xl0~ 4 to 
10~ 9 Mpc~ 3 , while the H 2 rotational line and the PAH 
11.25/im feature have luminosities from 10 38 to 10 42 and 
from 10 37 to 10 41 respectively. 

We do not find any significant difference in the line 
LFs of Seyfert l's and Seyfert 2's, which are indistin- 
guishable at the 2cr level. We also computed the line 
luminosity functions separately for the AGN l's and the 
AGN 2's, using our new classification, and find no differ- 
ences, probably because the shift of the 20 HBLR from 
the Seyfert 2's class into the new AGN 1 class has not 
a noticeable (statistical) effect on the new AGN 1 lumi- 
nosity functions, while it increases the uncertainties in 
the AGN 2 luminosity function. We do not present the 
luminosity functions with the new classification, because 
the AGN 1 sample size is 53 but the AGN 2 sample is 
only 21, making the statistics of these latter quite poor. 

Our sample includes all the galaxies selected at 12/xm 
that show evidence of Seyfert activity through optical 
spectroscopy. We are aware that some AGN activity can 
be detected in the mid-IR in optically unidentified AGN. 
Goulding & Alexander (2009), using a volume limited 
sample to D< 15Mpc, find that ~ 50% of the AGN they 
detect are not identified using optical spectroscopy. How- 
ever, these objects are typically starburst dominated sys- 
tems hosting modest AGN activity and have low lumi- 
nosities. The luminosity of the [NeV]14.32/xm emission 
line measured in these galaxies is about 10 37_39 erg/sec, 
which is too low to affect our luminosity function (see 
Fig. 16a). 

8. THE ACCRETION POWER IN THE LOCAL UNIVERSE 

Because the [NeV]14.3/im line is uniquely generated by 
the accretion process through the ionizing continuum of 
the central engine, its luminosity function can be used to 
measure the accretion power in the Local Universe. The 
first step is to determine the correlation between the lu- 
minosity in this line and the 19/im continuum luminosity 
that we measured in 59 of our Seyfert galaxies (see Table 
8). This relation is: 

LogL^Z = 0.9667 x LogL([NeV]U.32fj,m) + 4.3263 

and is shown in Fig. 21a, for each of the galaxies de- 
tected in the [NeV]14.32/im line and for which we have a 
measure of the AGN component. We can use this relation 
to convert the [NeV]14.32^m line luminosity function, of 
the whole set of galaxies measured in this line, into a 
19^m AGN luminosity function. We present this latter 
in Fig. 21b. The integration over all luminosities and 



over the volume defined by the average redshift of the 
12MSG of z=0.030 (RMS) gives a measure of the accre- 
tion power in the local universe at 19^im of 2.1 x 10 46 
erg s _1 . 

If we want to convert this monochromatic power at 
19/xm into the bolometric power due to accretion, we 
simply use the relations between the bolometric luminos- 
ity and the infrared luminosities published in Spinoglio 
et al. (1995). Taking the typical spectral index for 
Seyfert galaxies a(12 — 19/xm) = -1 and the 12^m vs 
bolometric luminosity correlation, we find a bolometric 
power of accretion of 8.0 x 10 46 erg s -1 , over the lo- 
cal volume out to z=0.03. We are well aware that the 
relation between 12/im and bolometric luminosity has 
been derived from large-beam (few arcminutes square) 
IRAS data, however both quantities are affected by pos- 
sible extended emission in an analogous way. Moreover 
all the Seyferts are dominated by AGN emission, which 
is almost point-like. In any case, a better estimate of 
the bolometric correction will be done using the nuclear 
fluxes of the 12MSG that we are planning to measure 
through the Spitzer IRAC and MIPS images (Spinoglio 
et al. 2009, in prep.). 

A more serious concern is that not all of the AGN 2's 
now emitting [NeV] are generating accretion power right 
now, i.e., in this decade, even though they evidently did 
in some previous centuries. This is offset to some de- 
gree by the missing population of "weak" AGN 2's which 
are not optically classified as Seyferts because their AGN 
component produces less than 50% of their infrared emis- 
sion currently (which we call "non-Seyfert-dominated"). 
The existence of these AGN 2's continues to complicate 
many discussions about " AGN" , since they can leave and 
enter this category depending on details of the detection 
observations. But since the 14.32^m luminosities of all 
these AGN 2's are well below the characteristic "knee" 
in the luminosity function, their inclusion or exclusion 
should not alter the luminosity integral substantially. 

We compare the derived accretion power with the 
power which originates in the starburst component of our 
Seyfert galaxies. We follow a similar procedure by using 
the [NcII]12.81/im vs starburst luminosity relation, as 
taken from the starburst percentage of the sample galax- 
ies from Table 8: 

LogLfH m = 0.9897 x Log L([N el I] 12. 81 fim) + 2.0198 

We derive a monochromatic power at 19/xm due to star 
formation and stellar evolution in Seyfert galaxies of 2.3 
x 10 45 erg s _1 , which corresponds to about 1/10 of the 
accretion power. This is quite reasonable, since the host 
galaxies of unambiguous Seyfert nuclei are only about 
one ninth of the 12 Micron Galaxy survey. In summary, 
because of their rarity, the set of galaxies we classify as 
AGN, namely those producing NeV line emission, gen- 
erate only a few percent of the total fusion power in the 
local universe. 

9. DISCUSSION AND CONCLUSIONS 

Our large sample of Seyfert galaxies, and the accompa- 
nying spectropolarimetric classifications have improved 
and extended the analysis of Paper I. In the IR diagnos- 
tic diagrams presented, we find that AGN 1-defined as 
having broad line emission of some kind-have high values 
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of ionization-sensitive line ratios, relative to the strength 
of star-formation tracers. In contrast, those galaxies we 
re-classify as non-Seyfert's have low ionization-sensitive 
line ratios and high PAH equivalent widths, and lie in the 
opposite regions of the diagnostic diagrams from AGN. 
The class of AGN 2's-those without broad lines-spread 
across both regions of the diagrams. 

The simplest and strongest version of the AGN Unified 
Model requires that any differences between Seyfcrt l's 
and Seyfert 2's should be due to the inclination angle of 
an obscuring structure (e.g. of a torus) with the line of 
sight to the observer, the Seyfert 2's being covered by the 
thick dust in the equator. However, in order to salvage 
this unification model from the results presented here, 
one has to explain why AGN 2's show a wider range 
of mid-IR properties, with similarities both to AGN 1 
and the non-Seyfert's. Perhaps this could derive from 
the particular geometry of the "torus-BLR-scattering re- 
gion" : the torus in fact can absorb the radiation not only 
from the BLR but also from the scattering region where 
the polarized emission could be produced, depending, for 
example, from the precise inclination angle, or the intrin- 
sic dumpiness of the torus. That is the basic suggestion 
of Heisler, Lumsden & Bailey (1997). However, to sup- 
press high-ionization emission lines we observe in the IR 
which are assumed to be isotropic- the absorbing struc- 
ture would also need to block out a significant fraction 
of the ionizing photons that would otherwise illuminate 
the NLR, even though this does not seem to happen in 
Seyferts with any detectable BLR. But once additional 
intrinsic differences other than just the viewing angle are 
admitted, the attractive simplicity of ATM unification is 
lost. 

On the contrary, however, there may be "genuine" 
type 2 Seyfert nuclei. These would be the half of our 
AGN2 category which have smaller proportions of AGN 
emission, and do not have analogs among the Seyfert 
Is. As we discussed, the remaining NLR emission we 
detect in the optical and infrared spectra of these true 
Sy2's may be the "fossil" remnants of a 10-100 parsec 
extent of gas that had been ionized by an AGN which 
was active in previous millcnia, but has been effectively 
"turned off" for the last several hundred years. Given 
their strong variability over all timescales, AGN which 
recently "turned off" must exist. Depending on the duty 
cycles of power generation through black hole accretion, 
they could account for up to half of our AGN2. This 
would then imply the existence of a comparable popu- 
lation of "recently turned on" AGN. These could stand 
out as having unusually strong broad line and near-IR 
continuum emission, relatively to their NLR strength. 
The most likely candidates for these "young" or more 
accurately, "rejuvenated" AGNl's are those Seyfert l's 
having relatively strong Fell emission and low values 
of NLR/BLR ratios such as [OIII]5007/H/3 (Boroson & 
Green 1992). Thus it could well turn out that long-term 
variability is at least as important as viewing angle in 
unifying the various observed classes of AGN. 

- The main results of this paper are as follows: We 
present the Spitzer IRS high-resolution spectra of 
almost 80% of the Seyfert galaxies of the 12/im 
galaxy sample, a total of 91 galaxies; 

- We adopted a spectropolarimetric classification, 



with the " AGN 1" class broadly defined to include 
both the Seyfert l's and the "hidden broad line Re- 
gion" Seyfert 2's, as detected through optical spec- 
tropolarimetry. All of our infrared diagnostics are 
consistent with Sy l's and HBLR Sy 2's belonging 
to the same single class. The AGN 2 class contains 
the remaining Seyfert 2 galaxies without polarized 
broad lines arc likely a mixture of weaker AGN 1, 
in which the BLR exists but has not yet been de- 
tected, and "true" AGN 2, galaxies that may not 
have been producing much hard ionizing radiation 
for the last several hundred years. Our AGN 1/2 
distinction, based solely on reddening-independent 
IR data, is supported by more of the data than the 
usual traditional classification scheme that divides 
the Seyfert l's and 2's based strictly on the de- 
tectability of BLR emission in direct optical spec- 
troscopy. It appears that the mid-IR observed 
properties characterize the AGN 1 as an homoge- 
neous class of objects. The mid-IR behaviour of 
AGN 2's, instead, shows objects similar to AGN 
l's and others more similar to non-Seyfert galax- 
ies. 

- Semi-analytic models based on the observed mid- 
IR spectra are effective in separating the AGN and 
starburst components in Seyfert galaxies. We find 
that for 31 AGN 1 the average AGN percentage 
contribution to the 19^m luminosity is 92.2±6.6%, 
while for 16 AGN 2 this percentage decreases to 
79.3±15.7% and for 9 non-Seyfert's is 67.6±17.2%. 
Although with large scatter, there is a clear trend 
of decreasing AGN strength from AGN1 to AGN2 
to non-Seyfert's. Ionization-sensitive line ratios 
can discriminate AGN l's from AGN 2's and non- 
Seyferts. The diagnostic diagram of [NeV]/[NeII] 
versus [OIV]/[NeII] provides a measure of the AGN 
strength in both AGN l's and AGN 2's. The AGN 
2 do not appear to be a homogeneous population 
with respect to starburst tracers, such as the EW 
PAH and the EW of [Nell]. Moreover, some AGN 
2 nuclei could either contain more ongoing star for- 
mation and/or be covered by more dust than any 
of the AGN l's. 

- The 0-0 S(l) H 2 rotational transition can be used 
to estimate the mass of the emitting regions. The 
average mass value is of the order of 10 s M©, which 
is in agreement with other estimates for Seyfert and 
starburst galaxies. 

- The density-sensitive line ratios (the [NcV] and 
[SIII] doublet ratios), in about 40% of the objects 
imply electron densities below the low density limit. 
Therefore they can not be used reliably to estimate 
the electron density of these regions. The simple 
interpretation that the line emission from the NLR 
and the HII regions can be heavily absorbed by 
dust even in these mid-IR lines (Dudik et al. 2007) 
is not confirmed from our data because the galax- 
ies of our sample do not show both line ratios af- 
fected at the same time in the same objects. In 
addition, the abnormally low ratios occur with the 
same frequency in all different types of emission 
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line galaxies. The other possibility is that some as- 
sumed atomic physics parameters need revision for 
both of these line ratios. 

- We derived for the first time the line luminosity 
functions for either classical Seyfert l's and 2's and 
also for AGN 1 and AGN 2. We do not find signif- 
icant differences between the various populations, 
in either the shape or the normalization of their 
line LFs. 

- The mid-infrared [NeV] lines are unambiguous 
tracers of the AGN NLR. We therefore use the 
[NeV] line luminosity function of all Seyfert galax- 
ies to estimate the accretion power in the local uni- 
verse within a volume out to z=0.03. We find that 
the power originating from accretion at 19^m is 
<~ 2 x 10 46 erg s _1 , about 4 times less than the 
bolomctric power. For comparison, the power re- 
lated to star formation and stellar evolution in the 
Seyfert galaxies population at 19/im is one tenth of 
that one from accretion. 

APPENDIX 

A. CLASSIFICATION: NOTES ON INDIVIDUAL OBJECTS 

NGC526A: has no detected broad lines, according to new data from Bennert et al (2006), however broad Ha 
wings were observed more than 20 years ago (Lawrence 2009, priv. com.). Because of the possible changes in 
BLR characteristics over timescaales of the order of years, we conservatively classify this object as a AGN1. 

MRK1034 NED1 = Akn80: this has not a Seyfert like optical spectrum (Osterbrock & Phillips 1977). We 
classify it as a non-Sy (IRS spectrum in paper I). 

MRK1034 NED2 = Akn81: only narrow lines, slightly wider than the instrumental resolution (Osterbrock 
& Phillips 1977). We classify it as a non-BLR as no polarimetric observations are available (IRS spectrum in 
paper I). 

ES0545-G13 = MCG-03-07-011 = MBG 02223-1922: detected broad Ha and H/3 ~ 2000 km/s (Coziol 
et al. 1993). We classify this object as an AGN1 (IRS spectrum in paper I). 

NGC1056: this object has been classified as an HII region galaxy through optical spectroscopy by Veilleux et 
al. (1995). 

NGC1097: detected broad Ha - 10000 km/s FWHM (Storchi-Bergmann, Baldwin & Wilson 1993). On the 
basis of this BLR evidence, we classify this object into the AGN1 class. 

NGC1194: only narrow lines were detected: Ha -230 km/s FWHM and [OIIi]5007A~400 km/s FWHM 
(Kirhakos & Steiner 1990). Because no polarization data are available, we classify it as a non-BLR object, but 
we are unable to include it in cither AGN1 or AGN2 classes. 

NGC1365: detected broad H/3 of 1896 km/s FWHM (Shulz et al. 1999). On the basis of the BLR evidence, 
we classify this object into the AGN1 class. Note: [NeV] is detected. 

ESO362-G018 = MCG-05-13-017: H/3 has a witdh of 5240 ± 500 km/s FWHM, therefore is a genuine AGN1 
(Bennert et al 2006). 

NGC2639: detected broad Ha of 3879 km/s FWZI (Keel 1983). On the basis of the BLR evidence, we classify 
this object into the AGN1 class. 

NGC3511: only narrow lines were detected: Ha ~ 135 km/s FWHM, this galaxy is classified as HII region 
galaxy (Kirhakos & Steiner 1990). We classify this object as a non-Sy. 

MCG+00-29-23: this object has been classified as an HII region galaxy through optical spectroscopy by 
Veilleux et al. (1995). 

NGC4602: detected only marginal evidence for an Ha broad component with FWZI— 10000 km/s (Kollatschny 
& Fricke 1985). We classify it as a non-BLR, as no polarimetric observations are available (IRS spectrum in 
paper I). 
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NGC4922: this object has been classified as intermediate between a LINER and an HII region galaxy through 
optical spectroscopy by Veilleux et al. (1995). 

NGC5005: this object is classified as a LINER galaxy (Goodrich & Keel 1986; Terashima, Ho & Ptak 2000). 

NGC5506: observed in spectropolarization: no broad Ha component was detected at a level of lxl0~15 cgs 
(Lumsden et al. 2004). We classify this object as an non-HBLR, and thus AGN2. 

NGC6810: this object is classified as a transition object between an HII region galaxy and an AGN, having an 
[OIII]A5007l width of 304 km/s FWHM and an Ha width of 263 km/s FWHM and the following line ratios: 
[Nil] /Ha = 0.62, [OIIIJ/H/3 =0.6, [SII]/Ha = 0.3 (Kirhakos & Stciner 1990). We conservatively classify it an a 
non-Sy. 

MRK897: Rodriguez, Malkan & Spinoglio (2009) confirm the classification of this galaxy as a Seyfet 2 of 
Durrct (1994). Because of the lack of polarization data, we classify it as a non-BLR (IRS spectrum in paper I). 

NGC7130 = IC5135: this object is classified as a LINER(Veilleux et al. 1995). 

NGC7314: detected in spectropolarization at Ha : F6563(broad)=5. 6xl0~15 cgs (Lumsden et al. 2004). We 
classify this object as an HBLR, and thus AGN1. 

NGC7496: the relative strength of the emission lines in the nucleus of this object is typical of normal photoion- 
ization found in HII regions (Pence & Blackman 1984) (IRS spectrum in paper I). 

NGC7590: the relative strength of the emission lines in the nucleus of this object is typical of normal photoion- 
ization found in HII regions (Storchi-Bergmann, Kinney, & Challis 1995) (IRS spectrum in paper I). 

CGCG381-051: this object has been classified as an HII region galaxy through optical spectroscopy by de 
Grijp et al. (1992). 

B. PECULIARITIES OF THE "20^M PEAKERS" 

We have considered the " 20/xm peakers" identified by Wu et al 2009 in section 2 and computed the average values 
of the discussed quantities in the classification Tab. 6. We find in our sample fifteen "20^m peakers": 11 Syl's, 3 
HBLR's and 1 AGN2. The major peculiarities of this group of objects are that: they are almost all compact at 19/ito, 
they have very faint PAH 11.25/ito emission, a flat average spectral index (60-25)/im, and low equivalent widths of the 
lines (see Tab. 6). 

In particular, among them 11 are compact sources (belonging to the extendedness class I), 1 belongs to class II, 3 do 
not have background subtracted spectra and therefore no extendedness measurement. 7 of them do not show any PAH 
feature at 11.25 /mi and the mean PAH EW 11.25 /im among the other 8 objects is -0.04±0.04/jto. The mean EW of 
the lines and the feature in Tab. 6 are consistent with the lowest values for the AGN1. On the other hand, the AGN 
contribution is among the highest values in the AGNl's, as it can be seen from the mean [NeV]14.32/im/[NeII]12.81/im 
ratio, the mean spectral index and the average modelled AGN contribution. 

C. MATHEMATICAL TREATMENT OF THE MODELS 

We present here the details of the semi-analytic models. To model the observed IR quantities, we have found an 
analytical expression for each of them as a function of Sft, which is defined as the ratio of the Starburst dust continuum 
at 19/xm to the AGN flux at the same wavelength. 5R varies from zero - emission totally from the AGN - to infinity - 
emission totally from the Starburst. 

pgal 

®=J^f (CD 

r 19fJ,m 

The extendedness is defined as the ratio between the flux measured at 19.5^im from the LH slit to that measured 
at 19/xm by the SH slit. If we define Ffg^ m and F^n as the fluxes of galaxy and AGN, respectively, through the SH 
slit, considering that at maximum the LH slit can intercept 4 times the Ffg and that F^n ^ oes Q0 ^ d e p en d on 
the slit aperture, as can be approximated to a point-like source, then: Flux(LH) = F$^n + 4 . Ffg^ m ; Flux(SH) = 

pAGN , pgal 
r 19fJ.m T r 19ftm- 
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We know that the PAH emission feature at 11.25/im originates in the host galaxy only, while its underlying continuum 
derives from both galaxy and AGN. Therefore we can predict its equivalent width as: 
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EWPAH = . 1 ' n - 25 :" (C3) 
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where Ai is a costant depending on the measured features flux and we adopt Ai = 1.8 for the PAH at 11.25/im. From 
the assumption that we can model the continua with a power law: 



v~ a F v = k (C4) 

where a is the spectral index, k is a constant. We adopt a=-lA for the AGN component and a— -2 for the galaxy, 
these values are in agreement with the values commonly used as spectral index for AGN and starburst and availble 
in literature (cf Edelson & Malkan (1986)). We derive the Ff^ 25Aim and the F^ G 2 ^ m as functions of Ff^ m and 
Fi9/4?n> respectively, and consequentely as functions of the free parameter SR. From equations C3 and C4, the PAH 
EW becomes: 

1.8- ( )(-2+i-4) . £ 
EWPAH = h ,, k J m (C5) 

\k/11.2bnm> JlTI 

The [Nell] EW has contributions both from the galaxy and the AGN: 

EW[NeII] = M-F^ m + B^U (ce) 

^12.8lMm + ^12.81/im 

Because the galaxy component is stronger than the AGN's in the [Nell] emission, we chose the constant A2 two 
orders of magnitude greater than B 2 , adopting A 2 =0.25 and B 2 =0.001. 

The procedure for modeling the line ratios is similar. We modelled the emission of each single line and then we 
simply took their ratio: 

[NeV]/[NeII] = M ' Hl }^ m + B * ' , (C7) 

A jpgal , R pAGN 

[OIV]/[NeII] = ^' F *?* m l* A ' F ™r , (C8) 

A 2 ■ ^12.81 M m + ^2 ' ^12. 8 l M m 

with A 3 = and B 3 = 0.01 and A 4 = and B 4 = 0.01, because we consider the [NcV]14.32/im and [OIV]25.89/zm 
lines to be produced exclusively by the AGN. 
To model the a( 60 _25)^m spectral index, we use its definition as: 

writing F ly (60/xm) and F„(25/im) in terms of F gal and F AGN , considering that F„(60/xm) is due mainly to the galaxy 
emission, it becomes: 

l0 9\ p9 al , p AGN ) 

r 25pm Tr 25(im //-| in \ 

"(60-25)^ - l oa ( v ™»™ \ ' ^ ' 

The treatment of the previous equations is the same of the EW PAH model's. All the factors A, and B^ have been 
found by fitting the data with the models and they depend on the instrument (spectrograph properties, slits areas, 
spectral extraction methods). These same models can be suitable for other observables and instruments, by scaling 
those factors. 



D. PAH 11.25/xM EW - LOW/HIGH RESOLUTION COMPARISON 

We plot in Figs. 22a and b the PAH 11.25/im EW measured by Wu et al. (2009) versus that one measured by us, to 
compare the results on the objects in common. In Fig. 22a we present the ratio of the values found by us to the values 
of Wu ct al. (2009) to better show the scatter of the two sets of data. In Fig. 22b we plot the PAH 11.25/im EW by 
Wu et al. (2009) versus ours and we computed their least square fits, calculated the mean slopes and their confidence 
intervals by using the bootstrap method (by 1000 resampling). The two sets of data are consistent with each other, as 
it is shown by the slope of the lines reproducing the least square fits of the data: considering all the galaxy types, the 
slope of the fitting line is 1.14±0.17; if we consider only the AGN1 and AGN2 the slope is 0.97±0.18. These two results 
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agree within the errors. However, the fitting line slope becomes steeper when we add the non-Sy's data, because for 
these we used a different method to measure the PAH emission feature at 11.25/im. Wu et al. (2009) remove the 
continuum underlying the range 10. 80-11. 80^m and integrate the PAH feature flux inside this window, in each of the 
sources. On the other hand, we measure the flux in the range 11. 15-11. 65/xm for all the sources, but we subtract the 
continuum underlying the whole range ll.2-13.0yum for those sources having in this range a broad multicomponent 
feature. This results in a lower continuum baseline measure and therefore in a higher PAH feature equivalent width. 

E. ADDITIONAL OBSERVATIONAL MEASUREMENTS 
Weak ionic fine structure lines 

Besides the ionic lines we presented in §4, the spectra of seven sources showed other fine structure lines, e.g. 
[ArV]13.10/zm, [FeII]25.99^m and 24.81/zm, [FeIII]22.95/im and 33.04/xra. The [FeII]25.99^m line is often blended 
with the [OIV]25.89/im line, because the IRS spectral resolution is not high enough to resolve these two lines in all 
the objects. These seven galaxies and their respective lines fluxes are reported in Tab. 12. 

There are unfortunately not enough good detections of these lines in this sample to base strong statistical tests 
on them. We have nevertheless compared our [FeII]25.99^m detections with those of Dale et al. (2009) in the line 
ratio diagram of [SiII]34.8/im/[SIII]33.5/xm vs [FeII]25.99/zm/[NeII]12.8/im. Although the small number of points, our 
measurements are in agreement with theirs, as the few AGN 2's are located in the region of the Seyfert's of Dale et 
al. (2009) and the non-Sy are located among their normal spiral galaxies. 

Full Width Half Maximum of the spectral lines 

The lines FWHM are due to the velocity field of the gaseous regions where they are produced. The velocity dispersion 
a of the [OIII] at 5007A, originated in the Narrow Line Regions, is related to the black hole mass (Nelson 2000; Greene 
& Ho 2005). It is interesting to extend the same relation to the mid-infrared lines of [NeV] at 14.32/im, [NcV] at 
24.31/um and the [OIV] at 25.89/^m, because they suffer less obscuration than the optical lines and are mostly originated 
by the AGN activity, though the [OIV] can be contaminated by strong starburst emission. 

Dasyra et al. (2008) claim to have measured resolved spectral widths of the IRS lines. The IRS resolving power 
for the high resolution modules is on average (500 +/- 50) km/sec along the whole spectrum. We find the FWHM 
of all the lines to be, within the errors, of the same value of the resolving element and to have a scatter of about 
the 15% of the mean value (the mean value and the standard deviation of the FWHM of all the fine structure lines 
are reported in Tab. 13). We cannot conclude that the lines are resolved, because we consider a line to be resolved if 
its FWHM is significantly larger than the resolution element. While Dasyra et al. (2008) consider the instrumental 
velocity dispersion (the one corresponding to the resolving power of the instrument) as the ratio between resolution 
element and 2.35, we compare the measured FWHM with the resolution element itself, to be more conservative. 

Moreover, both the typical AGN lines and those originated by starburst emission have nearly the same mean FWHM 
(see Tab. 13), which does not change with the Seyfert or galaxy type. The AGN lines, originated from the NLR, should 
suffer from more broadening than the lines emitted by the galactic regions. We conclude that the broadening of the 
lines is not intrinsic to the gas motion, but most probably is due to instrumental effects, otherwise we would measure 
different widths from lines originated from the different regions. 

In conclusion, we consider that the FWHM of the lines measured in the IRS spectra do not measure the dispersion 
velocity field in the NLR, and higher spectral resolution data would be necessary. 
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NOTES ON PROGRAMS: 1: Data without off-source measurements; 2: Data from Spitzer Archive; 3: Data from P50253 to be released; 4: see Appendix A for the source classification. 



TABLE 2 

Fine Structure Lines in our Sample 



Line fluxes (10 14 erg s 1 cm 2 ) in SH Line fluxes (10 14 erg s 1 cm 2 ) in LH 



NAME 
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TABLE 2 

Fine Structure Lines in our Sample 



NGC7130 


5.27±0.84 


79.3±0.93 


9.09±0.64 


29.4±0.77 


19.6±0.33 


54.2±0.88 


5.22±0.86 


19.7±0.84 


48.2±2.59 


93.9±4.90 


NGC7172 


5.87±0.61 


33.0±1.01 


10.2±0.67 


17.1±0.68 


11.9±1.00 




13.8±0.49 


45.4±0.48 


26.9±1.51 


59.3±2.42 


NCG7213 


2.32±0.61 


25.7±0.74 


<1.85 


12.0±0.80 


5.00±0.85 




<7.89 


<13.5 


6.97±0.88 


15.7±1.61 


3C445 


2.75±0.32 


2.31±0.24 


1.95±0.21 


6.23±0.01 




<3.54 


5.84±0.11 


22.6±0.92 


<3.16 


<3.69 


NGC7314 


15.9±0.53 


8.08±0.39 


16.9±0.52 


23.2±0.53 


9.97±0.71 




21.5±0.32 


67.0±0.41 


15.0±1.71 


14.2±1.76 


MCG-03-58-007 


2.75±0.21 


8.52±0.34 


6.63±0.33 


9.29±0.37 


5.20±0.60 


1.92±0.36 


3.91±0.30 


8.80±0.34 


<2.84 


11.7±0.95 


NGC7582 


21.3±1.43 


322±6.41 


38.8±1.54 


105±2.05 


87.3±1.99 




63.6±4.29 


262±5.54 


244±7.85 




NGC7674 


16.1±0.27 


20.1±0.3 


21.2±0.3 


35.3±0.32 




21.6±0.74 


16.5±0.46 


49.3±1.33 


14.4±1.22 


29.7±1.08 


CGCG381-051 


<0.61 


19.1±0.44 


<0.70 


1.35±0.32 


<0.65 




<1.18 


<1.26 


9.36±1.13 


12.5±1.72 



20 



Tommasin et al. 



TABLE 3 

Molecular Hydrogen lines and PAH emission feature at 11.25/iM 



Line fluxes (10 14 erg s 1 cm 2 ) 



NAME 


H 2 S(3) 


H 2 S(2) 


H 2 S(l) 


H 2 S(0) 


PAH 


EQ.W 




(9.67^ra) 


(12.28/um) 


(17.04/um) 


(28.22/im) 


(11.25/im) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


MRK335 




<0.16 


0.39±0.05 


<0.72 






MRK938 




2.72±0.38 


7.06±0.65 


<1.02 


292 


-0.687 


MRK348 




<1.05 


1.77±0.48 


<1.14 






NGC526A 




<1.65 


2.40±0.77 


<1.41 






IRAS01475-0740 




<1.20 


2.75±0.30 


<2.18 


30.3 


-0.071 


NGC1056 




<2.46 


3.58±0.48 


2.25±0.27 


214 


-1.196 


NGC1097 


<1.69 


6.95±0.55 


14.4±0.24 


4.95±0.94 


285 


-0.911 


NGC1142 


1.88±0.10 


1.09±0.10 


4.60±0.08 


2.76±0.21 


73.5 


-0.707 


MCG-02-08-039 




0.89±0.29 


1.21±0.38 


<0.64 


5.86 


-0.021 


NGC1194 




<0.83 


1.65±0.25 


<0.96 






NGC1241 




<0.83 


4.62±0.45 


<0.57 


44.7 


-0.453 


NGC1320 




<1.26 


3.26±0.37 


<1.14 


51.5 


-0.072 


NGC1365 




5.38±1.65 


21. Oil. 17 


37.59±14.87 


657 


-0.4 


NGC1386 




<2.75 


5.65±0.86 


<2.91 


84.1 


-0.106 


IRAS03450+0055 


<2.11 


1.20±0.29 


<2.12 


<2.18 






3C120 


<1.68 


<1.94 


<2.52 


<1.72 






MRK618 


<3.50 


<2.54 


2.41±0.48 


<4.38 


29.7 


-0.047 


F04385-0828 




1.90±0.42 


3.06±0.46 


<2.17 


50.5 


-0.060 


NGC1667 




1.49±0.31 


3.43±0.36 


3.66±0.48 


49.3 


-0.807 


ESO362-G018 




<0.94 


2.21±0.22 


1.59±0.20 


71.7 


-0.232 


ES0253-G3 


<1.08 


<1.02 


1.93±0.39 


<1.36 


46.7 


-0.092 


F05563-3820 




<1.22 


<0.92 


<2.88 






MRK79 


<0.99 


<2.19 


4.54±0.39 


<2.18 


13.4 


-0.020 


NGC2639 




<0.6 


0.96±0.34 


<3.33 


13.6 


-0.199 


MRK704 


<2.00 


<1.69 


<3.06 


<7.25 






MRK1239 




<1.13 


<1.29 


<0.85 


34.2 


-0.023 


3C234 


<1.48 


<1.75 


<0.57 


<1.51 






NGC3511 




2.11±0.17 


2.13±0.19 


6.53±0.46 


63.4 


-1.086 


MCG+00-29-023 


<2.26 


<1.55 


8.00±0.51 


<6.54 


259 


-0.702 


NGC3982 




<1.00 


2.65±0.36 


<1.16 


52.6 


-0.597 


NGC4051 




2.52±0.77 


8.79±0.54 


<7.37 


133 


-0.122 


MRK766 




<1.30 


1.78±0.35 


<1.45 


58.5 


-0.074 


NGC4388 




5.60±0.62 


14.5±1.27 


12.3±1.55 


103 


-0.141 


NGC4593 




<0.98 


1.71±0.40 


<4.51 


44.5 


-0.051 


NGC4922 


<1.69 


1.98±0.35 


7.81±0.47 


<2.70 


68.2 


-0.145 


NGC4941 




<1.73 


4.61±0.47 


<1.35 






NCG5005 




11.96±0.84 


26.0±0.62 


8.35±0.85 


243 


-1.240 


MCG-03-34-064 




6.96±0.54 


6.98±0.80 


<2.75 


71.2 


-0.044 


MCG-06-30-015 




<1.69 


2.91 ±0.52 


<5.94 


25.5 


-0.035 


NGC5135 




<1.14 


4.97±0.25 


8.94±0.81 


124 


-0.554 


NGC5256 




1.63±0.19 


4.68±0.24 


<3.33 


41.0 


-0.595 


IC4329A 




<2.06 


<3.19 


<6.83 






NCG5347 




<2.12 


6.49±0.91 


2.95±0.73 


33.2 


-0.060 


NGC5506 




1.74±0.33 


3.75±0.24 


<3.40 


34.8 


-0.053 


NCG5548 




<1.70 


<2.45 


<0.25 


34.9 


-0.085 


NGC5929 




<1.03 


4.10±0.27 


<1.37 


22.6 


-0.495 


MCG-02-40-004 




<1.37 


6.80±0.38 


1.79±0.48 


91.2 


-0.120 


FSC15480-0344 


<0.89 


<0.80 


0.99±0.32 


<0.98 


49.2 


-0.142 


NGC6810 




3.81±0.59 


14.2±1.06 


<3.50 


91.2 


-0.486 


NGC6860 




<1.56 


2.69±0.24 


1.18±0.24 


21.7 


-0.056 


MRK509 


<2.25 


<1.24 


<2.62 


<1.79 


53.9 


-0.093 


IC5063 




6.64±0.97 


9.43±1.03 


8.65±1.42 






NGC7130 




<2.64 


8.01±0.72 


<2.94 


216 


-0.455 


NGC7172 




<2.13 


6.60±0.76 


4.23±0.67 


145 


-0.410 


NCG7213 




<1.73 


2.45 ±0.76 


<4.31 


62.2 


-0.116 


3C445 


<0.91 


<0.61 


<0.72 


<2.70 






NGC7314 




<1.01 


1.99±0.46 


<2.00 






MCG-03-58-007 




<0.77 


3.07±0.44 


<1.11 


47.1 


-0.081 


NGC7582 




9.06±2.02 


25.6±2.70 


<14 


920 


-0.490 


NGC7674 




1.64±0.25 


7.18±0.22 


<1.54 


84.6 


-0.091 


CGCG381-051 


<1.49 


<1.89 


3.80±0.43 


1.93±0.41 


43.8 


-0.203 



Columns (1) to (4) give the line fluxes of the H2 rotational lines in units of 10 14 erg s 1 cm 2 . Column (5) and (6) give the flux and the 
equivalent width of the PAH 11.25 fim emission feature, respectively. 
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TABLE 4 

H2 Estimated Temperatures and Masses 



NAME 


T(2-3l 


TCl-21 


T(0-1) 


H2 Mass 




(K) 


(TO 


(TO 


no 8 m^o 


MRK938 




304 


>404 


<1.47 


MKN348 




<379 


>184 


<0.51 


NGC526A 




<414 


190 


1.05 


IRAS01475-0740 




<322 


>171 


<1.32 


NGC1056 




<413 


186 


0.12 


NGC1097 


<195 


339 


206 


0.27 


NGC1142 


363 


249 


189 


4.65 


MCG-02-08-039 




432 


>198 


<1.23 


NGC1194 




<346 


>191 


<0.37 


NGC1241 




<226 


>462 


<11.93 


NGC1320 




<304 


>237 


<0.21 


NGC1365 




257 


134 


1.93 


NGC1386 




<340 


>265 


<0.03 


IRAS03450+0055 


366 


369 


157 


3.97 


MKR618 




<563 


>134 


<9.72 


F04385-0828 




389 


>178 


<0.98 


NGC1667 




322 


155 


1.57 


ESO362-G018 




<318 


177 


0.48 


ES0253-G3 




<355 


>178 


<5.02 


MKR79 




<339 


206 


2.35 


NGC2639 




<390 


>114 


<0.68 


NGC3511 




535 


117 


0.15 


MCG+00-29-023 




<232 


>169 


<0.22 


NGC3982 




<301 


>214 


<0.04 


NGC4051 




304 


>168 


<0.08 


MKN766 




<430 


>170 


<0.46 


NGC4388 




304 


167 


1.17 


NGC4593 




<372 


>122 


<0.60 


NGC4922 


<277 


256 


>138 


<12.12 


NGC4941 




<301 


>258 


<1.02 


NGC5005 




331 


>246 


<0.20 


MCG-03-34-064 




<538 


224 


1.71 


MCG-06-30-015 




<374 


>129 


<0.63 


NGC5256 




291 


>178 


<5.02 


IC4329A 




408 


128 


2.93 


NGC5347 




<283 


211 


0.39 


NGC5929 




<255 


>242 


<0.23 


MCG-02-40-004 




235 


>273 


<3.08 


FSC15480-0344 




<459 


>162 


<1.63 


NGC6810 




261 


>283 


<1.63 


NGC6860 




<373 


214 


0.58 


IC5063 




420 


>163 


<2.08 


NGC7130 




<284 


>231 


<1.79 


NGC7172 




<281 


184 


0.37 


NGC7213 




<421 


>135 


<0.52 


NGC7314 




<348 


>158 


<0.08 


MCG-03-58-007 




<255 


>233 


<2.65 


NGC7582 




<293 


>196 


<0.80 


NCG7674 




245 


>307 


<3.09 


CGCG381-051 




<368 


192 


3.80 
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TABLE 5 

Photometric fluxes and extendedness factors 



NAME 


Fl9 Mm (SH) 

(Jy) 


F 


l9.5um(LH) 

(Jy) 


R 


Class 


ivi rirvyoo 





742±0 


.006 


n 
u 


.885±0. 


.009 


1 19 









,589±0 


.007 


n 
u 


.647±0. 


.003 


1 in 

J — LU 


j 


MPP^9fi A 





363±0. 


.017 


n 
U 


361±0. 


.003 


u.yyo 


i 


1IX.tj.oVJ ±'-t I o-u I *±\J 





,462±0 


.003 


n 
u 


.507±0. 


.004 




I 


IN Uv^lUOD 





130±0 


.006 


n 
U 


.209±0. 


.001 


±.D± 


TTT 
ill 


NPP1 1 AO 





096±0. 


.001 


n 
u 


.184±0. 


.015 


1 92 


III 


mpp 09 ns n^o 

iv± v_< vj - u z - u o - u o y 





356±0 


.005 


n 
u 


379±0. 


003 


1 nfi 


J 


MP Pi 1 





336±0 


.007 


n 
U 


.427±0. 


004 


1 97 


TT 
11 


MPPl 941 





,100±0. 


.004 


n 
U 


.144±0. 


.002 


1 AA 


TT 
11 


IMP PI ^90 





698±0. 


.005 


n 
u 


.743±0. 


.005 


l.UD 


1 

1 


IN *aVylOOU 





793±0 


.015 


n 
u 


990±0 


.006 


1 23 


JJ 


1 ix/\ o U o^O U ~r u u 





,473±0 


.006 


n 
U 


.463±0 


.003 


u.yo 


T 
1 


QPl Of| 





,579±0 


.024 


n 
U 


.577±0 


.006 


u.yyo 


1 

1 


MRKfil & 

1V1IXJA.O-LO 





,568±0 


.012 


n 
u 


,577±0 


.026 


1 02 









910±0 


.007 


1 

1 


.008±0 


.006 


i.ii 


T 
1 


TMPP1 fifi7 

IN LtV-/10u 1 





053±0 


.004 


n 


.154±0 


.009 


9 Q1 
z.yi 


III 


t^qoiro Pms 





339±0 


.005 


n 
U 


361±0 


.002 


l.UD 


T 
1 


Tr , Q09^Q pq 





,565±0 


.003 


n 
U 


667±0 


.002 


1 1 8 
1. lo 


TT 
11 


9^Q 





,944±0 


.001 


n 
u 


.971±0 


.002 


1 07 
l.U i 


J 







,237±0 


.005 


n 
u 


,258±0 


.005 


1 HQ 
i.uy 




MPP*3Q89 





.116±0 


.003 


n 
U 


.172±0 


.003 


1 A% 


TT 
11 


MT? 'Wlf^P. 
IVlIX-iY 1 00 





,873±0 


.002 


U. 


999±0 


.003 


1 1 Q 

l.lo 


TT 
11 


MPPA^SR 


1 


.225±0 


Oil 


1 


.445±0 


.008 


1 1 S 
l.lo 


TT 
11 


MPPAQzH 

IN vj V_> i ±y i ± 1 





209±0 


.020 


n 
u 


.247±0 


.006 


1 lo 

J LO 


JJ 


mpp ^nn^ 

IN L'VjDUUO 





139±0 


.015 


n 
U 


.240±0 


.005 


1 79 


TTT 
111 


MPP D 1 } n^Q fifizL 


1 


990±0 


Oil 


J 


.992±0 


.005 


1 on 

l.UU 


I 


mpp ofi Qfi ni ^ 

IVI V_< Vjr- U O- O U- U 1 





681±0 


.019 


n 
U 


679±0 


.010 


u.yyo 


T 
1 


MPP ^^AS 





,425±0 


.006 


n 
U 


.502±0 


.005 


117 
1.1 1 


TT 
11 


MPP 09 ACi C\f\A 





,576±0 


.006 


n 
u 


690±0 


.005 


1 2 


JJ 


rqpi c;A8fl fi^zLzt 





,455±0 


.006 


n 
U 


496±0 


.002 


1 HQ 

i.uy 


1 

1 


MPPfiRi n 


1 


,514±0 


.005 


J 


910±0 


.004 


1 26 


JJ 


NGC6860 





,280±0 


.003 





.317±0 


.004 


1.13 


I 


MRK509 





,557±0 


.007 





.631±0 


.004 


1.13 


I 


IC5063 


2 


,344±0 


.008 


2 


.480±0 


Oil 


1.06 


I 


NGC7130 





,814±0 


.009 


1 


.094±0 


.004 


1.34 


II 


NGC7172 





256±0 


.014 





338±0 


.003 


1.32 


11 


NCG7213 





,425±0 


.020 





541±0 


.003 


1.27 


II 


NGC7314 





233±0 


.008 





.248±0 


.004 


1.07 


I 


MCG-03-58-007 





,651±0 


.005 





,677±0 


.005 


1.04 


I 


NGC7582 


2 


337±0 


.042 


3 


.290±0 


.015 


1.40 


11 


CGCG381-051 





336±0 


.008 





.343±0 


Oil 


1.02 


I 
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TABLE 6 
CLASSIFICATION TABLE 



Syl 



HBLR 



AGN1 



AGN2 



non-Sy 



"20/im peakers" 



Extendedncss 
EW PAH 11.25/xm 

[NeV] 
pVeTTJ 

EW[NeII]12.81/im 
EW[OIV]25.89Mm 
EW[H 2 ]17.07Mm 
%AGN at 19/im 

"(60-25)^™ 



1.17±0.29 
-0.15±0.20 

0.60±0.51 

-0.028±0.027 

-0.122±0.127 

-0.009±0.008 

92% ± 6% 

-0.84±0.62 



1.17±0.22 
-0.16±0.19 
0.73±0.53 
-0.028±0.026 
-0.082±0.069 
-0.008±0.006 
92% ± 8% 
-0.86±0.81 



1.17±0.26 
-0.15±0.19 
0.65±0.52 
-0.028±0.026 
-0.116±0.134 
-0.009±0.008 
92% ± 7% 
-0.86±0.74 



1.39±0.45 
-0.35±0.27 
0.66±0.84 
-0.081±0.057 
-0.093±0.054 
-0.046±0.069 
79% ± 16% 
-1.65±0.93 



1.54±0.58 
-0.81±0.42 

0.10±0.03 
-0.116±0.029 
-0.033±0.035 
-0.058±0.062 
69% ± 16% 
-2.19±0.65 



1.08±0.08 
-0.04±0.04 

1.03±0.87 
-0.014±0.015 
-0.075±0.047 
-0.008±0.003 
96.8% ± 1.5% 

-0.44±0.65 



Note: the mean values have been computed without considering the upper limits for all quantities and the zero values for the PAH EW. 
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TABLE 7 
Equivalent widths of key lines 



NAME [NeII]12.81/xm [NeV]14.32/im [NeIII]15.56/im H 2 17.02/^m [OIV]25.89/m?i 





H 


(H 


W 


W 


(H 


MRK335 


-0.002 


-0.005 


-0.010 


-0.005 




H TT) TV noo 

MKK9oo 


-0.106 




-0.025 


-0.025 




MRK348 


-0.030 


-0.008 


-0.037 


-0.004 


-0.063 


JNGG526A 


-0.017 


-0.022 


-0.040 


-0.013 


-0.165 


1KA&01475-0740 


-0.028 


-0.006 


-0.023 


-0.008 


-0.023 


NGC1056 


-0.103 




-0.093 


-0.025 


-0.007 


T\T /"< /~1 -1 -1 J o 

NGG1142 


-0.141 


-0.016 


-0.063 


-0.056 


-0.032 


MGG-02-08-039 


-0.011 


-0.022 


-0.030 


-0.006 


-0.103 


T\T /"< /~1 1 -1 f\ A 

JNGC1194 


-0.007 


-0.009 


-0.018 


-0.005 


-0.060 


TVT /"I /"I 1 O ,1 1 


-0.113 


-0.015 


-0.076 


-0.045 


-0.066 


JNGC1320 


-0.013 


-0.019 


-0.022 


-0.006 


-0.066 


t\t/~i t~\ 1 o o i - ■ 

NGC1386 


-0.019 


-0.040 


-0.045 


-0.007 


-0.154 


tt> A OAO/icn i nn cr cr 

!KAb03450+0055 


-0.001 




-0.008 




-0.010 


3C120 


-0.017 


-0.038 


-0.065 




-0.461 


TV T T> TV I o 

MKK618 


-0.028 


-0.007 


-0.009 


-0.004 


-0.031 


F04385-0828 


-0.015 


-0.001 


-0.009 


-0.004 


-0.013 


JNGC1667 


-0.183 


-0.023 


-0.195 


-0.078 


-0.065 


EbO362-G018 


-0.043 


-0.012 


-0.026 


-0.007 


-0.044 


ES0253-G3 


-0.028 


-0.014 


-0.040 


-0.004 


-0.058 


T~y a r r n o OOOn 

1 05563-3820 


-0.005 


-0.004 


-0.009 






MRK1239 


-0.007 


-0.002 


-0.010 




-0.037 


3C234 




-0.010 


-0.017 




-0.099 


JNGG3982 


-0.100 


-0.039 


-0.073 


-0.028 


-0.050 


TV TDT/^r^ 

MKK766 


-0.032 


-0.028 


-0.031 


-0.002 


-0.071 


TXT/"^ 1 -i o o o 

JNGG4388 


-0.088 


-0.054 


-0.116 


-0.014 


-0.296 


T\T/~1 /~1 , 1 n A 1 

JNGL-4941 


-0.075 


-0.050 


-0.143 


-0.028 


-0.225 


NGG5005 


-0.136 




-0.091 


0.178 


-0.044 


TV /I no on c\c a 

MGG-03-39-0o4 


-0.034 


-0.039 


-0.070 


-0.004 


0.1088 


MGG-Oo-30-015 


-0.010 


-0.007 


-0.009 


-0.007 


-0.076 


JNGG5135 


-0.176 


-0.032 


-0.087 


-0.021 




T /~i ■! o nn A 

1G4329A 


-0.014 


-0.014 


-0.031 






JNGG5347 


-0.006 


-0.004 


-0.008 






JNGG5506 


-0.034 


-0.029 


-0.072 


-0.007 




NCG5548 


-0.020 


-0.013 


-0.017 




-0.076 


TV if /""i /~i n n /in nn *i 

MGG-02-40-004 


-0.022 


-0.010 


-0.015 


-0.013 


-0.030 


r bG15480-0344 


-0.017 


-0.018 


-0.028 


-0.002 


-0.120 


NGC6810 


-0.086 




-0.014 


-0.013 


-0.005 


NGC6860 


-0.015 


-0.010 


-0.023 


-0.010 


-0.087 


MRK509 


-0.027 


-0.006 


-0.032 




-0.096 


IC5063 


-0.012 


-0.015 


-0.033 


-0.005 


-0.075 


NGC7130 


-0.146 


-0.016 


-0.052 


-0.012 


-0.020 


NGC7172 


-0.092 


-0.026 


-0.046 


-0.019 


-0.163 


NCG7213 


-0.053 




-0.031 


-0.007 




NGC7314 


-0.039 


-0.082 


-0.106 


-0.007 


-0.467 


MCG-03-58-007 


-0.013 


-0.013 


-0.017 


-0.005 


-0.024 


NGC7582 


-0.126 


-0.022 


-0.060 


-0.016 


-0.088 


CGCG381-051 


-0.087 




-0.006 







* MRK335, F05563-3820, NGC5135, IC4329A, NGC5347 and NGC5506 do not have calculated [OIV] equivalent width, because even if they have 
background subtracted spectra, the continuum in the (19-35)/J.m range results not well calibrated. 
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TABLE 8 

AGN and Starburst Contribution at 19/ai 



NAME 


AGN % 


Starburst % 


error 


Seyfert 1 


ESO012-G021 


85.4 


14.6 


4.5 


ESO545-G013 


78.4 


21.6 


13.5 


NGC0931 


97.7 


2.3 


0.7 


MKR618 


95.6 


4.4 


3.4 


ESO362-G018 


85.9 


14.1 


4.3 


MRK0006 


98.2 


1.8 


0.9 


MRK0009 


97.6 


2.4 


0.7 


MKR79 


98 


2 


1.1 


MKN1239 


97.7 


2.3 


1.6 


NGC3516 


95.6 


4.4 


1.1 


MKN766 


92.8 


7.2 


3.9 


NGC4748 


81.4 


18.6 


15.2 


MRK0817 


92.1 


7.9 


5.0 


IRASF15091-2107 


92.5 


7.5 


3.4 


ESO141-G055 


97 


3 


1.7 


NGC6860 


93.9 


6.1 


3.4 


MKN509 


90.9 


9.1 


4.6 


HBLR 


ESO541-IG012 


95.9 


4.1 


1.9 


NGC0424 


96.9 


3.1 


2.2 


NGC0513 


72.5 


27.5 


3.2 


IRAS01475-0740 


93.8 


6.2 


2.9 


NGC1125 


79.1 


20.9 


5 


MCG-02-08-039 


97 


3 


0.5 


F04385-0828 


95.8 


4.2 


2.3 


NCG4388 


91 


9 


3.1 


NCG5347 


92.9 


7.1 


4.1 


TOLOL01238-364 


84.4 


15.6 


3.7 


MCG-03-39-064 


96.9 


3.1 


2.0 


MCG-02-40-004 


94.6 


5.4 


0.8 


FSC15480-0344 


94.5 


5.5 


0.5 


IRASF22017+0319 


95.6 


4.4 


1.3 


MCG-03-58-007 


95.7 


4.3 


1.7 


AGN 2 


MRK938 


52.9 


47.1 


12.6 


IRAS00198-7926 


96.8 


3.2 


0.8 


NGC1142 


44.8 


55.2 


6.2 


NGC1241 


69 


31 


7.9 


NGC1320 


95.4 


4.6 


1.4 


NGC1667 


71.7 


28.3 


11.9 


ESO033-G002 


96.8 


3.2 


1.9 


ES0253-G3 


89.6 


10.4 


2 


NGC3660 


75.5 


24.5 


7.3 


NGC3982 


80.6 


19.4 


7.7 


NGC4501 


70.2 


29.8 


7.6 


NGC4968 


90.1 


9.9 


2 


NGC6890 


89.5 


10.5 


1.3 


NGC7582 


68.8 


31.2 


11.3 


NCG7172 


85.3 


14.7 


3.5 


non-BLR 


MRK1034 NED02 


63.8 


36.2 


17.5 


NGC4602 


68 


32 


11.9 


non-Sy 


MRK1034 NED01 


67.0 


33.0 


21.0 


NGC1056 


68 


32 


12.1 


NGC5005 


43.1 


56.9 


8.5 


MCG-03-34-063 


41.6 


58.4 


17.7 


NGC6810 


80.5 


19.5 


8.8 


MRK0897 


53.8 


46.2 


16.1 


NGC7130 


77.2 


22.8 


9.8 


NGC7496 


81.2 


18.8 


9.3 


NGC7590 


75.6 


24.4 


10.3 


CGCG381-051 


87.5 


12.5 


14 
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TABLE 9 
S-K Statistic for AGN families 



Diagram Syl + HBLR Syl + non-HBLR HBLR + non-HBLR 



F(H 2 ) vs F(PAH) 
L(H 2 ) vs L(PAH) 



[Nelll] 
pVeTT] 
[NeV] 
[Nell] 
[OIV] 



Nell] 
NeV] 



EW PAH vs 
EW PAH vs 
EW PAH vs 
EW PAH vs jgjjjj 
EW PAH vs EW[NeII] 
cxt vs EW[NeII] 
ext vs EW[NcIII] 
cxt vs EW PAH 

[Nelll] 
jNeTlT VS 01 
[NeV] 

WeTT] VS a 



757o 


1 007 

187o 


77% 


79% 


88% 


20% 


100% 


15% 


39% 


54% 


23% 


54% 


40% 


8% 


41% 


98% 


9% 


20% 


78% 


0.1% 


0.4% 


60% 


19% 


10% 


69% 


17% 


19% 


84% 


14% 


33% 


97% 


12% 


12% 


100% 


6% 


5% 



TABLE 10 
Line Luminosity Functions of Seyfert 1 



Luminosity (erg s x ) 






Luminosity Function(10 


6 Mpc 3 mag 


X ) 








L 


[SIV] 


[Nell] 


[NeV] 


[Ncill] 


[SIII] 


[NeV] 


[OIV] 


[SIII] 


[Sill] 




(10.51/im) 


(12.81/um) 


(14.32^m) 


(15.56/iin) 


(18.71/mi) 


(24.32^tm) 


(25.89/im) 


(33.48/j.m) 


(34.82/im) 


38.5 
38.9 


13.43±11.05 




16.28±11.41 














39.3 


44.29±30.75 


38.14±30.44 


41.44±30.62 


38.14±30.44 


38.14±30.44 


54.43±32.50 








39.7 


44.50±25.28 


40.23±24.95 


46.74±25.35 


34.97±20.21 


71.09±28.68 


28.14±17.51 


43.43±30.62 


168.69±95.09 


143.33±93.51 


40.1 


10.17±6.91 


27.81±13.77 


10.17±6.91 


25.20±18.86 


2.83±2.26 


31.81±19.85 


16.72±11.53 


9.86±7.11 


25.31±17.26 


40.5 


4.07±2.05 


7.11±4.06 


1.61±0.75 


11.05±6.94 


6.05±2.70 


1.22±0.67 


44.82±25.28 


21.30±12.02 


30.26±13.94 


40.9 


2.15±1.11 


17.61±2.62 


4.68±2.20 


3.45±1.95 


2.03±1.11 


4.94±2.26 


10.58±6.93 


6.95±2.82 


3.70±2.01 


41.3 


0.36±0.22 


1.95±1.11 


0.27±0.22 


1.02±0.46 


0.83±0.39 


0.75±0.38 


3.42±1.95 


1.10±0.48 


4.17±2.07 


41.7 




0.16±0.09 




1.62±1.06 




0.03±0.01 


1.76±1.09 


0.26±0.01 


1.05±0.47 


42.1 














0.12±0.07 




0.16±0.11 



Line Luminosity Functions of Seyfert 2 



38.5 


114.28±81.44 




32.52±29.47 






81.76±75.92 


81.76±75.92 






v . 

c+ 


38.9 


86.32±64.58 




166.88±99.73 


2.76±2.00 




81.13±64.54 






32.52±29.47 


N 
CD 


39.3 


64.12±48.31 


100.48±88.84 


72.75±49.44 


128.04±82.38 


42.74±31.89 


109.23±57.71 


120.74±71.16 


1.70±1.70 


1.66±1.34 


t~i 

CO 


39.7 


99.55±88.65 


71.17±49.72 


108.10±89.08 


182.06±109.84 


342.08±142.49 


115.61±89.25 


20.78±12.50 


244.56±120.38 


112.97±89.73 


40.1 


28.34±9.46 


219.44±104.26 


19.11±7.77 


23.77±13.23 


18.98±7.03 


21.65±8.42 


75.87±49.19 


137.51±82.39 


171.40±99.82 


CD 
& 


40.5 


4.60±2.16 


10.16±6.42 


19.07±7.95 


17.41±7.76 


20.93±7.91 


8.89±4.65 


3.96±1.84 


29.19±10.96 


39.86±14.89 




40.9 


2.40±0.88 


15.14±6.29 


2.18±0.88 


8.91±3.41 


3.12±1.12 


2.00±0.82 


26.22±9.35 


11.37±4.96 


10.80±5.97 




41.3 


0.90±0.74 


2.48±1.02 


1.11±0.77 


2.36±0.86 


0.89±0.42 


1.53±0.88 


1.01±0.44 


1.35±0.59 


3.07±1.53 


C 
^ 


41.7 




1.68±0.64 


0.14±0.08 




0.12±0.08 




0.52±0.44 


1.02±0.51 


1.28±0.59 


CO 


42.1 














0.29±0.17 


0.07±0.03 


0.81±0.48 


CD 
'< 
CD 


Luminosity(erg s _1 ) 


Seyfert 1 


Seyfert 2 












~t galaxi 


L 


H 2 S(l) 


PAH 


H 2 S(l) 


PAH 












cd 

CO 



38 


5 




20.87±11.50 




137.56±76.19 


38 


9 


56.82±35.52 


3.04±2.84 


197.81±111.08 


20.16±12.41 


39 


3 


42.54±20.89 


0.38±0.20 


66.75±58.40 


3.66±1.80 


39 


7 


12.41±7.36 




55.29±19.55 


2.56±1.41 


40 


1 


5.90±2.93 


0.31±0.25 


15.10±5.01 


1.30±0.54 


10 


5 


2.34±1.12 


0.09±0.08 


1.44±0.70 




40 


9 


0.72±0.32 




2.55±0.80 




41 


3 










41 


7 










42 


1 











to 
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TABLE 11 

New Classification, PAH flux and EQ.W. and H2 Masses for Sources 

in Paper I 



NAME Old class. New class. PAH at 11.25^m H 2 Mass 

flux EQ.W. (10 8 M ) 

(lO^ 14 erg s^cm,- 2 ) (/im) 



IRAS00198-7926 


Sv2 

yjy a 


non-HBLR 


31.2 


-0. 


066 


<19.36 


ESO012-G021 






88.9 


-0 


417 


<1.91 


IRAS00521-7054 


Sv2 


HBLR 








<12.27 


ESO541-IG012 


Sv2 

uy ^ 


HBLR 


12.5 


-0 


031 


<27.80 


NGC0424 


Sv2 


HBLR 


19.2 


-0 


001 


<0.42 


NGC051 3 


Sv2 


HBLR 


74.7 


-0 


588 


3.69 


MRK1034 NED01 (1) 


Svl 


Tl OTl-Sv 


87.1 


-1, 


105 


5.36 


MRK1034 NED02 (1) 


Svl 


non-BLR 


187 


-1. 


051 


<11.86 


ESO545-G013 (1) 


Svl 




21.9 


-0. 


262 


2.37 


NGC0931 


Svl 




39.5 


-0 


043 




NGC1125 


Sv2 


HBLR 


64.2 


-0. 


465 


<0.33 


ESO033-G002 


Sv2 


non-HBLR 


10.8 


-0 


024 


0.66 


MRK0006 


Svl 




20.0 


-0. 


041 


<30.29 


MRK0009 


Svl 




15.2 


-0. 


039 


<2.20 


NGC3516 


Svl 

oy ± 




25.8 


-0 


037 


<0.22 


NGC3660 


Sy2 


non-HBLR 


32.4 


-0. 


584 


0.94 


NGC4501 


Sy2 


non-HBLR 


36.1 


-0 


652 


0.41 


NGC4602 (1) 


Syl 


non-BLR 


28.2 


-0. 


689 


0.38 


TOLOL01238-364 


Sy2 


HBLR 


111 


-0. 


131 


<0.44 


NGC4748 


Syl 




20.7 


-0. 


.300 


0.38 


NGC4968 


Sy2 


non-HBLR 


88.1 


-0 


133 


<0.30 


MCG-03-34-063 


Syl 


non-Sy 


16.7 


-1. 


176 


<1.24 


MRK0817 


Syl 




24.5 


-0 


043 


<2.73 


IRASF15091-2107 


Syl 




25.7 


-0. 


081 


<9.36 


ESO141-G055 


Syl 




19.2 


-0. 


056 




NGC6890 


Sy2 


non-HBLR 


48.8 


-0. 


155 


0.20 


MRK0897 


Sy2 


non-Sy 


117 


-1. 


194 


1.66 


IRASF22017+0319 


Sy2 


HBLR 


9.56 


-0. 


021 


<6.75 


NGC7496 


Syl 


non-Sy 


153 


-0. 


.419 


<0.12 


NGC7590 


Syl 


non-Sy 


56.1 


-1. 


166 


0.12 



Sec Appendix A for the source classification. 
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TABLE 12 
Rare ionic fine structure lines 



Line fluxes (10 14 erg s 1 cm 2 ) 
NAME [ArV] [Felll] [Fell] ' [Fell] [Felll] 

(13.10/im) (22.93^™) (24.52/um) (25.99^tm) (33.04^tm) 



NGC1056 




2.83±0.04 


NGC1142 




2.51±0.03 


NGC1194 




0.85±0.01 


NGC1241 




2.41±0.01 


ES0253-G3 


0.02786 




NGC5135 




6.80±0.61 


NGC6810 




11.12±0.05 ■■• 9.69±0.12 1.23±0.08 


NGC7130 




11.81±0.10 
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TABLE 13 

Average Full Width Half Maximum of the Lines 



Dispersion Velocity (km s 1 ) 





b WriM 


St. Dev. 


[SJV]10.51/im 


685 


122 


[NeII]12.8lixm 


604 


151 


[NeV]14.32[tm 


625 


123 


[NeIII]15.56fj,m 


638 


119 


[SIII]18.71fim 


673 


185 


[NeV]2A.32fj,m 


601 


163 


[OIV]25.89^m 


552 


89 


[SIII]33A8tim 


580 


119 


[SiII]3A.82^m 


708 


215 


H 2 S(l) 17.04,um 


587 


149 
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MRK335 



o.i 




Mrest) (Mm) A(rest) {/im) \(rest) [//m) 




IRASO 1475-0740 



S(3) S(2) 



,f 



t 1# 



/ 1 



jr 



H a S{1) 17.03(j.m 



□ 



[SHI] 33.48/im 



v [Nell] [Nelll] 

[SrV] [NeV] [SI 



20 

X(rest) (^m) 



NGC1056 



il, s(i) tr.oaiim - 




[OIV] 35!B9/jm 1 




. [Fell] _ 








V 


17V. 



0^ 




[Nell] 


[Nelll] 


[NeV] 


[sm] I 


[NeV] 


[sin] 




[OIV] 



20 

X(rest) [/im) 



Fig. 1. — Spitzer IRS SH and LH specra of the observed Seyfert galaxies. Wavelengths have been shifted to the galaxies rest frames. For 
the objects with no off-source observation the SH spectrum is shown, because slightly affected from the background emission (lsimeq 10 
%), together with the > 3cr detected lines in separated boxes. Figures 1.11.61 are available in the online version of the Journal. 
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NGC1097 
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Fig. 2. — H2 excitation diagrams. For each measured line, the natural logarithm of the level population normalized to its statistical 
weight is plotted against the upper level energy (in temperature units). For each pair of adjacent transitions the connecting line is shown, 
whose inverse value represents the gas temperature: the dashed line connects the S(0) and S(l) detections, the dotted line the S(l) and 
S(2), the solid line the S(2) and S(3). Upper limits have been used to obtain limiting slopes and hence limiting temperatures and masses 
(see text). Figures 2.12.43 are available in the online version of the Journal. 
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Fig. 3. — a: H2 17.04 (im line flux versus PAH 11.25|tm integrated flux, b: H2 17.04 fim line luminosity versus PAH 11.25/im luminosity. 




Fig. 4. — a: The PAH 11.25/im equivalent width versus [NeV]14.3/^m/[NeII]12.8/^m line ratio. The black line shows the behaviour of the 
analytical model for this diagram, b: The PAH 11.25^tm equivalent width versus [OIV]25.9/jm/[NeII]12.8/xm line ratio. The black line as 
in Fig. 4a. 
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Fig. 5. — a: The PAH 11.25/im equivalent width versus [NeV]14.3/mi/[SiII]34.8/jm line ratio, b: [NeIII]15.5/im/[NeII]12.8/im line ratio 
versus the PAH 11.25/xm equivalent width. 
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Fig. 6. — [NeV] 14.3/im/[NeV]24.3/im line ratio versus the [SIII]18.7/^m/[SIII]33.5/im line ratio. Symbols as in the previous figures, except 
for the open crossed symbols indicating the ratio of objects with small aperture [SIII] 18.71^tm measures and the filled symbols directly 
above them showing the aperture corrected [SIII] line ratio. The dashed lines show the low density limits (see text). The solid lines at the 
top and at the right give the corresponding electron densities. 



50 



Tommasin et al. 



AGNl - SYl 
non-BLR 
AGNl - HBLR 
AGN 2 

non-SY, HII, SB gal 
non-SY; LINER 



G 



6. 



log(n) = 3.25 
* • • ^ 

low density limit ^^^^ 



en 



G 



(a) 



AGNl - SYl 
non-BLR 
AGNl - HBLR 
AGN 2 

non-SY, HII, SB gal 
non-SY; LINER 



A 



6 



low density limit 



TL? Mm m 

- .. . £* • 

• • ■ ■ !■ 



og(n)=3.0 



(bX 



[SIV[10.51/im]/[SiII]34.B2^m] 



[NeV[14.32/im]/[NeII]12.81/im] 



Fig. 7.— a: [SIV]10.51/xm/[SiII]34.82/im line ratio versus the [SIII]18.7/xm/[SIII]33.5/Ltm line ratio. The [SIII]18.7/mi flux has been 
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Fig. 8.— a: [NeV]14.3/im/[NeII]12.8^m line ratio versus the [OIV]25.9/mi/[NeII]12.8/mi line ratio. The black line shows the behaviour 
of the analytical model for this diagram, b: [NeII]12.8^tm equivalent width versus the PAH 11.25/^m equivalent width. We note that, for 
graphical reasons, in this diagram and in the following ones in which the equivalent widths, covering a large range, the logarithm of the 
inverse of the actual equivalent width is plotted. The black line shows the behaviour of the analytical model for this diagram. 
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Fig. 9. — a: PAH 11.25/im equivalent widths versus extendness parameter. The black line shows the behaviour of the analytical model 
for this diagram, b: [NcII]12.8/tm equivalent widths versus extendness parameter. The black line shows the behaviour of the analytical 
model for this diagram. 
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Fig. 12. — a: [NeV]14.3/im/[NeII]12.8/xm line ratio versus the 60-25/xm spectral index. The black line shows the behaviour of the 
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Fig. 13. — The histogram shows the percentage of AGN contribution. 
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Fig. 14. — Average high resolution spectra for our classes of galaxies, compared with the mean high resolution spectrum of starburst 
galaxies. 
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Fig. 15. — a: Line luminosity function of the [SIV]10.51^m, the red squares represent the Seyfert l's and the blue triangles the Seyfert 
2's. b: Line luminosity function of the [NeII]12.81/im, symbols are as in the previous figure. 
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Fig. 18. — a: Line luminosity function of the [OIV]25.89/im, symbols are as in the previous figure, b: Line luminosity function of the 
[SIII]33.48/im, symbols are as in the previous figure. 
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Fig. 19. — Line luminosity function of the [SiII]34.82mum, symbols are as in the previous figure. 
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Fig. 20. — a: Line luminosity function of the H2 17.04/Um, symbols are as in the previous figure, b: Line luminosity function of the PAH 
11.25/jm, symbols are as in the previous figure. 
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Fig. 22.— a: The ratio of the PAH 11.25/im EW calculated by Wu et al. (2009) to ours as a function of the values by Wu et al. (2009). 
b: PAH 11.25/xm EW calculated by Wu et al. (2009) versus ours. We computed a mean square fit among the AGN1 + AGN2 class (values 
are fitted by the red line) and all the galaxy types (black line). For the slope values and confidence intervals, see the text. 
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